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2011 BAAS El X E SUPE HB TEE 


шен 
(EAE, Жй 100084) 


2011 JE [R^ RC НЕ НИ ЖЖ. КАНЕ — КИ ЕЖЕН Тау tS ЛИНЕ 
# T 1992 0175, MA 1994 fF иа аи aj A [R T. qp ts IRE REA Se ETA A 7, 2002 f 
Жн CH LS MOMS “у ХЕШ. ЕДЕН, Made ATL 
Ж” ЖМИ Е, ЕВИ ТАГО A OR a BS tš F FRI 
FARM. 20° FEAT, ТАЯТА В 3140 вт: WE, SHR 
BORGES 1251 97, BEST RI AT RS ИИ — №: BENRA 19490 А (® HEF 58000 £ 
Л), BRE EHECR ECS} ЈИ tg tF EJ Ct 16% R 24% ERK. MOREE RAR 
ARIA AGRE RAS AL, WÑ K. WA ZRELA, ІП ALR ОКА ЕУ 
MER, пи TR АЫ ЛЕНО. ШЕ EROS НА ИЕ CH OAT RL, 

2011 # Ft H 4 [RC 2F ^E BC: НЕ (ID 98 20 Мі Й) T 2011 F 9 1914512 AAT. 
BL 33 4 $6 / HAK (0446 fat RUN IAE) DL ДИ m R TRE 1251 Bii Et 19490 ^^ 
WA ОС ЖЕНДІ 16008 BA, FHL 3248 BA) HEKE 6 KA pin ЧЕН, EKS 
RAREN, 55 2010 F ff) 1196 ТЕЖКО 17317 BA (Xo ЖН 14119 BA. БИН 3198 HA) AE 
H, НИК 4.6%, B BOB 12.5% (Ж ЕНД RUP pf K 13.4% 18.0%). HERARET ЈЕ. 
JONG SEA AAT RLS РЈ OED Ж, JEEE T Matlab BIK. 

JUI fp W %% eh d I Ii PC 6 RMT OR ARE (^E HAR. PR. W 
№, ИМИ пин, WIB SORDO CRURA E RESELL A E EE 
AP PAL FE Ad S & (GE. мижеияжа. Жаян ШК “SKEER 
ЖИШШ ЖЖ”, СИ ИЙТ Ж). AE IF FI e XO WC 1821 0 (A, B, C, D IB 4} H 
№ 1037, 463, 204, 1768). fü SPE MH 9%. А BETI 1 HEIR О ORY. MA 
Жем Тану енен, укинтищех, 9 epo T НА. 
TOS. BLS AEG KO ELLA (SEHR. ни, БРОНЕ, вм 
PATE LC OPE AED Зимы, МЕНА 1536, ӨЗЕНІ 8%, 
SEERA YT. кожна-жзад, 0936104004. НИН 936510. = 
8 %221 a ЗЕН EES Sr f) 15% #16.5%. ЖЕБЕ КН ЖЖЖ 
ТЕХ. BEKE. ща. ИНТЕ AES. TEMERE. RR 
Ж. Matlab IRR ЖЕНУ B= ERA RR. KIM. ак, CRAMMER 
тш. иж. же. мм, жж. ЖЖ. Ш. Ев. Ж. жп. шж. ГЖ. як. М 
INF ASB EAT BR THER. 

AERE ВТ ЗЕЕ ИНН, УМЕНИЕ И: RM Т 
Ж. FAAS. ЗЕРЕН РНЕ НЕ ји LAR. ЗЕН А e Ee II BOE RI 
ЖЕШИНЕ ЖОО SER IU RU, BE SE. RW. HRS. писин, кални 


2 E туа = = KW ast 


ШААНИ; EXT EAUR- B. ИЕ, AREE = ЖИ УЕ виж 

X. Өр РИНЕ. еще, женя, елін. НКМ ШЕ ВЕ (НЕ AES 
мәннен AH CFE): RRM EER. ЕЗИНЕ. вина E 
пипай, ЖТИ ДӨ NATE. 

ута упавжамиет, (еи НЕ. ЖЖЖИ ЕИ "RU 
хета, CERT НА ЖЕШ {ЕМ AA (ева, ГЕШ. M EHS) А 
ЕНЕМ À (PET) век. HE RETIRE. MEKE, FETERE, 
жағал RAF. AEE. IUE RMT A, Випава" . 

АКА "oto ЕО" dili BET K ME HEH, ВМ “Ж 
AREP fa ЈЕ ОМ, E E L ЕС isa FUTURI ЕАУ A 
ваши, си “ФАШ КИТЕ $ MT JF” Hik KW ИН, DIE 
“KARR тиш” нин Жы МНН їн. ЕЕ ЖЖ. HTE Z. i 
ROLE A ЗЕЕ Е НОВУЮ, ЕТК ЕЕ A EL ване НАШЕ, 
ПОСЕТЕ aa e pK as ASI AR GW СУ ТОРЕ ПАРУ REAN 
# (Email: zhijiecai@163.com), ЭКЕЙ f ЗЕ (Subject) PHEW “КМ”, ТЕМЕ 
Фра, ке. ME, ШЖ. рим, ЕВЕ, US ЖМИ ХИНЕН 
WES MA Af 

AT SP ХИН, ERE LE S, РСА ЖИА Т Ie US. 
HES (EAS FR (НИВУ OER тания” Indie X). 
ИФА UCU. 

жн LIVEN ктш, № ТА: TE ON НСС КИЕ, X 
жарава БН ЖР АСА 

ВИЉА ЈЕ ВЕНЕ TIE T MURS. ШЫ ЯЛЕ: http://mcm.edu.cn. 


АШ WOSERTRESBRISHAHI 


BE M ti 6 A 2 ОВУ A CL AS Ди, A DISCS И: БЕЗИ ТЕЧ 
НН. APR ЧЕ ES IS жена атулар ана W u YE FY FF НЕ 
Rot, яла F Жозе ек, Нав Лижие. 

HME RIS, ЖКГА. THE. Ше. ЕЕ АВИВ, % 
Желі ни. 20. .... SRK, ЖИЙ АЖАЛ КА. 

ЖАР Ж RC + SMB FIG KATUR ft. e db. PUS С ТУ D LIS LC 
ПОН DOM, Н арг УНЕ AR REE (0-10 BOK ED) ЖИН. 64%. Ж 
HI GPS RRR AMEE. BUH PTK НА. заета ЧЕЖЕ ГЕ HLS CRO 
WORD. AT, BIRD 2 RR ДЫ VERE Е ЗА ARERR. КРЕ 
ОСЕ Шия то 

үкімін TREE РАНИ ЛЕ. На ЈЕЗЕРА НЕР HEN, MEIH T 8 Е: 
Bock B MRE, ЮЗИ T ВЕ ЕЕЕ ФИЛ ЖЕЛИШ. 

ЗЕЕ ЖИРЕ BRERA FIER - 

Е Е 
ва: 

2) HUET. ҮЙӨ ЯФ РДЕ TEKEL 


СІҢ шел: юпдинетнхтъменияя з 


3) MESAHE, BIE TEN, WER R 

Жалшалзал лалы А, XE EXIGERE MET В ЗЕЙ ЛЕН. жик. 
йв? пите УРА Г 

На НИЕ, Ж http:// mem.edu.cn Fa. 


BE СИНА РЕНЕ 5ME 


“WARS. йк WT. SER OUT ALA. ти, ELI 
Ж. MIME. H TERAM, ЖЕНЫ ESE ЗН НИИ 
W W Р Зр т. О ® 6 IT E УТИ ITA та 
+, ЈАДАР СА TOR POE РОТИРА эк 6 0 e 
MW. USE A DE IAE BE TaN] KEEL. 

ARCH ал WIB 5 Y E XISL, SY A ТІ ГА FII 09648; 

1) Mt 1 h fP MS Н ОБС А ВЖЕ IR CREDIT (0 20 I Fe M 
аллы HUC RON MIRE 2. Жз DERM FE ЕТІН, мат 
Matin dH Ра HR E PET порта зна fi 26 ЗЕ (ЕН Ж Эл бо km/h) 810599. 
юл. 


ИРЕК ЧИ, WIWER 20 EF УЖЕ, Xe i EC 13 fete 
LIGUE DOR MW. IER 1 6083 EHI AAU. Meri ЕН 
fr mH Ж. 

БГ TERRE IM AY 6709 T ff ЖН Ay i И кеи, КЕ 
BRM 2F SFE, iE A BEN P f NETO. 

2) HEIN (EINE A, B,C, D, E, F) ЛАЙ, BER E W EU ER ЗР бнт 
% апа ҮЛЕ LG P fr ib W yK (UL PE) m шшш ACM, W 
Wiwa k. 

ЛААНЫ СР (08 32 416 80 ERE ТЕХНИ, fex ar OE, qul 
мерт. Dy T GERM MER. шығара SE Ag бо Зе Vi TT IG 
жж. 

WHEL АРЕНЫ НН SF RM АЕН. 

М2: SSE 55 V i ИЙИН. ЖЧК tcs t T fd). 

JE. МИНИ, VA http://mcm.edu.cn FAR. 


CH PURARLHES eR 


HER Uo RG, до AHL лала ИЕЛИК Я 
Di EGRE A НЕЯ, АЕ ORR ACM IG g. 

RN MIT AS HEE OR ktp “ная” “лер” ARAB, ET 
HEIL T W 8 MOS ен (20%) RUH AE 2 ЕНЕ), Петля ee He 
Bl (990) MAEDA. RATKARNS RARE RIS. жп, RTZEN AE 
Л GRAE) MT 5 Ht Во T VEZ (НИМ). тре SUB AK — ERER FI T EE 
ИЖ, АИ PRENE (RIE), META Еф де He IER (A 
KARES) EAGT A IES. ЗЕРНА ТЕКЛИ. wx 


4 гв ж т та аве 


може, ERP SARL, PAR ИЕ. РА PEFR, 
АЛИ ЕЕК а, HANER, HEREN IH. 

жеш инт. {ЕЕ ТА ЊЕ РЕТТЕН НЛ ЖЖ. ка XC БМ 
жшж. E30 F RRO RIGS. CR RL АСРИ TOR Ж. 
тәне. RI LOOK DE, АОН ВЫ ЛЯ REF LA ITER 
BERRY. 

жшн наа ЖЕ РИИ. HER TAN MEEN HER ME e 
WEA, CHAAR DMEM Lhe. жа ORI k ~E 
EBORE E UR. CAR ЛЕДЕНИ AL 09 E eo A Be EL ЗОЖ 
EDU TIE QI RARO. ЖЕ ЖЕ AE (EG Н FF HUF OSE 58.5%. кешік, ша 
CLARE REA Ж ЈРЗ бр NE ВОЗЕ МР ЫТ ДН aio rais КИН, 
КҮЛЛЕ MARO. JA ORI Ж @ ЛУН 
ета IE 

тш A ГЕН) T ШК. Bie 2 Ë 2009 4F Ж ЖУКИ BURL (UIT. 
ПЕРВАН. ПОВЕЗ ЈЕЛИ ФИД НИКЕ. иивожоквия, Wu FH: 

HM, RKP SF ТОН ЛЕ УМ НИКА EL, emma. HBA 
тері, КОЩА 2011 4:55 2035 Ш ЖӨНІН ТЕЖЕ 

FM. НУМИМЕЗ УЕЙ 2009 4F Uie PUT T БАН. Ime 
заживее ТИ ОНЕКИ Б ни, АЕ oe ELT. (12000 ERAMA 30 H ~ ЕЯ 
жаа, — ИЩИ (55 07,60 3,65%). ИН КӨ ОНА. 

IS (ПА ЧЖИ. 2000 ед 30 % емен ОРИЧ (55 Y , 60. 
3,059), HABWE- (WO 75 HE. HERRE e (0M CH GL. ЖИИ 
ARNT SUD HERE SB DVR, Mt r 000 2 taa k p ИИ JE @ 2 a k Y Ñi. 

мн, шаятын E A M FORE. НЕ ОО pt Y. dici ar ЖШ 
DIM UILLAM 

d: НЕЕ, ИА http://mcm.educn FAR. 


DH RMR 


ква (BL MIRI RC) Ети EEN, [I PILAM MRE 
SADNESS ESP а P (ЖМ). ЖАТ. BEEF ARMM Ат, ша 
жн rin W. VERUM PHEIN HATS СЛЕ AT нида (В). 

ПОСЛЕ. BER AOS KA IS, ш. ЗЖЗАЖЕЗЖИЯ, 3.53.91 3.5K 
ин. KRE. жале клада, KERK, оо ARAL Ей, usc 
Fa ey Ж. 

ATH EPOR. pinin 2, ж сайын SL OR AE RINK 
JOCURI. HURL AUR. Qi RUE TAMER AH “MAM 
А" TE. 

ASSURE RATE ЕД: 

1) жена AIR, Rell RAR а: 

2) жану Ж. AURORA, JURE: 

З) EREBE, SKM ЖЕ 20.5 ЖИВЕ, викат Ни 


mia HOR, гола те врстан 5 


4) REME SE ROB ID RR НВС, PY LBB. MKB жз CA RIE fr 
TT Ж-із5 RTL, НОЈЕ 1-7 Ж-із.5 жене 

5) АТЫН, ERLE 30 A ER. 

WHEL EEF, азжини, НЕТ ВНЕ ATE RI, а 
HRRK. 


Ку жама 
HAKAR) ЖАК №) жа Вл (Ж) 
з 65 20 во 
7 135 8 ю 
14 oc 5 89 
жо жиш. 

FAH] зза | 35-89 | 444 | 4549 | 554 | 5559 | 664 | 6560 
жи 43 59 39 4 27 28 м 21 
XAQO| т-та | тэ-т» | вва | 85-89 | 994 | 95-99 |10104|105109 
жа 24 24 20 25 2 23 21 18 
KA (K) | 11-114 [1155-119 [12-124 [125-129 [13-134 |13.5-13.9 | 14-144 [14.5-14.9 
кк 3 23 22 59 18 25 35 29 
ЖА. (Ж) [15-154 [15.5-159[16-164 16-169 17-174 | 17.5-199 | 18-184 185 189 
жк | зо | a | за | a | 45 | | so | u 
KA (K) [19-194 [195-199 |20-20.1 [205-209 | 21-214 215-219 |224 [225-229 
ж 52 63 49 35 27 16 12 2 
ЖА (ж) | 23-234 [2315-239 [24-44 | 265-249 |25-25.4 [255-259 
жа 0 6 0 0 0 1 
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ВЫШ. Инн т О Ев сЕ, ИНН 

f. US RUN. бернн И. У H ТРОА 
atas sms. 


2 Ван AE О atakay E 
Изикежнана ния, шия, О “ә” HERK, "en b TE, 


"YMK, AXAR., “ж” HARRER. RAE AA A AA E AR 
WB OR ROC LAR А, BIRA CIL OR ML 


Hl HARE 


4 r жете пан 


Е О: 
жине. oO P ERAS, комата. DH ТТІ? 
ии. JU F+ 

RRRA Руси), WARN (= 1,2, n), рп НАК. INT KK 
мшу), ШАНА (та), ЈЕНЕ 

A. 0 Sr = Ont, 
(ту = x a 1 
(у) Бај ба ке а) 


ие 
пети? а 
РМ Р, БЕЛИ. 
3j — PARE RIO CAL ВИСО ОА ЛА Kriging МИР, у ЕВЕ РУАНО ЕЕ, ЖІЛІК 
Mut Rt AR u y y iñ RII НА 5 42 ЖТ К КИК МИЛЕ ПОПИС T E ERU 
вини. 
WE si P(x, y) КЕМ 


ћу) e m+ Y Хеб), (3) 
= 
Jt P RET P, Ж МИЙИК r, (2) КАНН, тА, (i = 1,2, n) МАКИ, ЖА 


(4) 


АНИ: 


huy) = М, #=1,2,- 
ШЕЙНИ m MA, (4 = 1,2, -.- ,n) ОВНА ВИД 


(5) 


AA =b, (6) 
же 
Фи) фы) ° |n) 1 № [1 
(ги) Фта) с Фа) 1 м м 
Жы 5 АА, О A: Ры О СЕ: 
(rni) Фа) 52 Фм) 1 А hn 
1 № +. 100 т 1 
жа 


тө = У —z; + (и, (8) 
РЫ Р, SZ IE. MILII FER. GERN т А, (i = 1,2, n) ИА, ЖАН 
Ваша (3) тегін OE A(z, у) 

Fl PERS ЕН ИЕ СУРЕ л ЖИЫ а. 


Е RESON 45 


ралады 


У. 
TRESS RITTER IR. AP TREE! а. 

оза WARREN 
HET AR. ХИ Muller НЕ. ЈЕ Muller F 1949 ВОН В UCD HO R SS e 
ШҮК. ARAB Т PR ЛУ ЕВА ЖАЙ RADA, TO LE HERE ANED 
RR, ЯЖ, KIA Г g ç W fifty RECA. ПОН ОГЫ МИИ 
SATA, HEE AA ia Wait RE p u. Е 

Tyo = loga ga @ 
СОИ, иес С СЕСЕ ТУГЕ 
M15 
PION ВИИ Ls MESE BIO XA HIR У. 


Жа: акене. 


закана ялкелкек зала. 
Tom <0 о ІЗ 
OS lee <1 1 НАЯ 
1S lew < 2 TAR 
4 251. <3 3 FARR BAR 
BS leo <4 4 ала 
i AS <5 5 | miam 
Пе 25 в | мана 
$2 ASM 


MT We A Dy Mk Me РОНЕ, ЖН ОР ТО О ЗО 
VR ek rV t HJ ра НАРНИЯ, ik AE He MERI Re КИ DU Е 
ЕРТТЕП ЕТЕРІ ТТЫ ЖЕЛИ. За Dues Las 21928 I TERRI 
юни Kt, тне X 


и P= (10) 
oett ок ШЖ 261259; 
Аз: ABFROTRRRA DIR 
даға Р ЕГІЗ ЕТІ 
1 P <07 Es Яв 
п O7<P<10 жал ай» 
| ш 10<Р<20 r= жн 
у 20<Р<30 EL 221 
у P230 - жз» Früh 
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а RRAHIN t 


ал MRED 

ута, SBA CSM Е ЈЕ ПОЛА. ШЖ ФДУ 
ЖЕЙ ЖЮ) ИЖЕ, BA CANTO M AE ФЕЈ, БИЙИГИН TAA, ана атана алар 
Miri] DL ДЕН ее E MR. 

НОСЕ 
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е — 


p= u) 
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42 мжам 

HX FRM Юон F. ОАЕ Я РАТУ аа 
ЖИ MIE. 

WER SPAT AI S ME FAY REO ЖП, EN RRR. ЕЛИ 
ANDRA ИО E TOC PTE. таттан нт, во R ЖЖ, А, 

BARI НАПОИЛКИЛЖИЮАХ, fj d AUR: 

жаз iE WEZ ИИИ ЖЕ сит) (4, j = 1, тв 

PRS HORDETERCK IO MOS GER A 

жиа ARS IE RZ HIC t: 

PRS рза, ИВО ЮЖ. 

ISDA CRT ELI" ЖИ ЖЖ! ЖЕ. Ha), aa, зат Abi, Das Бо № 


жей рїї. ie 
к= (aisi as) Y= (nib 2 (12) 
№ 
Ma, = B(z — Ez) (z — Бг), Му = Ely = Ey) (y – Бу), аз) 
Mj, = May = E(z Ез) (y = Еу). (04) 


А, a, МЖК ма мама Mys EWER, ЗС At RRE AIT OUR. Ж 
Ханум ARRON 


1 
sen- La 


TUEN: Sim on МН, SUR MAAK ЖОЛЫ, Шу” = 


Био Гати уа 


ЭТАН ЕШ. REDS ROE. ке ЕЗ алан Hh 
ЗЕРИ ЕВ DRE ALE по ОУ ABLES. 

ddu = ulz, y t) ЖЕ ЖН? T Р(х, y) ЬЕ Te RAE, ЖИЛ GE HR, TA 
Sp RI 


ди š 
+ Угато. (16) 


ян вал REIMER aT 


ФОНЕ RB ti KS, ОЖ SA E (Katika, Вет 
ЭЙЕ q E ЖИЕ ВИИ h ӨРШІЛ, В а ос Vh. 
ЖЖ, ВВЕЛ q TS RTE RUR IR FIS GER, а x u, БОК, T4 


а= -kuVh, an 


ЖЧ! SORTS MOT EA ME PK K ЕЕ, к (> 0) HEARN. 
¥ (17) ЖҚА (16) Ж. A 


ди 
Sp BVA: Vu - kuV?h = 0. as) 

КИН Т M LO & WN FLETES TCR OREA шщ (i 1,2, ,n), ETT E 
tT unm) = ш, i=1,2 n (19) 


HABET TE: WS HAASE (18)- (19), BUSTA АЕ ERES S ВСИ [10]. ХА 
Mox gr. TI BLM ER ЖЫ. 

IH FHA M Et MOI. EERIE ЖЕ КИЙ, МН fep pte, WU 
$ kt, Шу (18) roto (очна я, (ае) 


ЕЈ = Уп. Vu uh = 0. (20) 
(MLL, cl EJ = 1. 
ANIR (20) ҮННЕН E 2 
УМ (1) 
Mitre k t 
t=Tix=xr, (02) 


ПАНОНИЈЕ ЖЯ (Т, хт) УНЕ x(t; xr). РАНЕНО, ЕЛИ ЛОИ и та: 


а 
i ram. (03) 
иий 
t=Tiu=ur, (24) 
PRESET A MOS 
т 
u(x,t) = urexp (— І Vhx(rixr))dr). (25) 
‚№ (21) #1 (25) НЕ ИН, AEF (век, 358 RETI СЕНИ) 
xml = x" + мж"), (26) 
ш" = exp ( зА умом), (27) 


2 = 


Jot v Dy T M WS ДЕТЕ кво миене. 


48 r w m... = LE 


PERL HEH ч x° (CHE ао, RI ih (26) Ж (27) SE ELAR ЖЕ ее < T) HS AE HEAL 
8 хо REG E QUERER u(x, tix). gh. ЗНАЕТ ЗЕЕ МИНИ P, EGENT < 
туяр E кожа 


Bi = (9140, 5) > U), (28) 
SEE REY Ju E fe fi) ER F ВСС. 


6 FR PAN 


EDERAL P, PUERO SEI S ИИМК VER. ERE TEX e 3c 
Wk, ЖСК C Matlab @) ИЕЛІ XA. HIME RMT EE. 

ПОВЕСТИ И NCBI БЕ Т KAKER, G K£ RUK 
ЖАРАР» МУРА НИЖНИ. EAH T Matlab ‘Pf griddata iñ t k НЕИН, Vei Mt 
OET A PHIL M APAA ELEM. 2} KER FE A FU ЛА EE f RII 
жи, алашы, EARP, PIR A ASR RE A 

НЕШ КЖ К. пшн OEE FI ЖСК A В PL 
№, {БЕЛИМ rok: ТИЛИНИ, жив Ar HAP RI. ШЕЛІ TACITO IRL +í 
ЖАЙ, TAB HEFA FI, AJ PETIA KHL DG А. 

РВИ ТАСТ ВВ Зонд да Т Hof MTA HAR ee hE. ГАА о АЕ 
POR FAD LEON A. GARR, UH АНЕ OE ERNIE 470900028, (151 
ТОНО. 

HOA OA Ба, ELA TE FIFRA RC ЖЫН. ИЗЕЛ BUE 
ГОТ ЗАЛДА ней. Н Hr HH W @ алаштан 
ИН, пазен T EBORE AR RN. mitapi, ICOM ЕЕН И ду +20), ЈЕ 
p gy N НИКО ОКНЕ ВИ, ЛЛ V CE HIDE HH RM 
жик eb s. 

шола аа 


доста 
бе (29) 


ДЕТТІ НИЛИ р, o AAA A RRO PSS ЕТЕ ЗН 


1 = log, max (ë = 1, 1). (30). 


3) МОНЕ ЛИНЕАРНИ HEHE Hie DUCK NEGO а 
AMAL n. МИРАЖИ КА 


Picus 


Жара REZM ZE. 

4) ТЕМЕЉИ НИНЕ ПАША ЖУ KEN ST инт. акапежият, EX 
EJ PH CAS RUIN. BERERE, RED ПОРЕ ЗНАНА ВА, XR 
КЕЙТ MT A НО RURUM, G ВЕЛИ РАНИ ЛЕ Я. COH SCIO НЕ АНИ 
BOR. ПИРЕ О. DUE, БЕ SMR 
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ЖИЕ REE, EOE T FARAH, OR PS МЕНЕ БИ ЯНЕ Ж SII a 
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problems. We also give some comments on the students’ papers. 
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Abstract: According to the scattered features of the given data, à triangle spine interpolation method i 
constructed based on the Delauney triangular subdivision technique. Ву utilising this method, the spatial dis- 
tribution for each heavy metal is then obtained. Based on the fact that the background value corresponda to an | 
interval for each heavy metal, the factor assessment model and the comprehensive evaluation method are estab- 
lished. "The assesment result for cach subregion is also provided. Besides, by using the multivariate clustering | 
and correlation methode, the relationships between different heavy metals is obtained, and the main reason 
underlying the urban heavy metal pollution is evaluated. To recognize the source location of the heavy mental 
pollution, three models are constructed respectively. The first one is established based on the triangular spline 
interpolation technique, the second based on the characteristic of convection-dispersion, and the third based on 
the distribution characteristics of convection or dispersion. By combining these three models, the method for 
recognizing the source location is proposed and Ия effectiveness substantiated by experiments, 

Keywords: triangle spline interpolation clustering analysis; correlation analysis; convection-dispersion equa 
tioni normal distribution fitting; evenly distributed fitting 
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Analysis of Heavy Metal Pollution in Urban Topsoil 


AN Chuar-kai!, HUANG Xiang-yi?, LIAN Chen? 
Advisor: DONG Zi-jing? 
(1- School of Electronic Engineering and Computer Science, Peking University, Beijing 100871; 
2- School of Mathematical Sciences, Peking University, Beijing 100871) 


Abstract: During the rapid pace of the Chinese industrialization in recent years, the environmental side effect, 
especially the soil pollution, becomes more and more serious. Accumulating many data of the вой pollution, 
researchers pay increasing attention on unveiling the relationship between human activities and the environment 
by using statistical and machine learning tools. This paper evaluates the pollution distribution and analyzes 
the pollution level of eight major heavy metals in а certain urban area by utilizing Kriging interpolation and 
principal component analysis in part I. In the hierarchic analytic way, in part 2 the paper finds out the main 
reasons (factors) underlying soil pollution and substantiates the relative weights of these factors. 

Keywords: Kriging interpolation; principal component analysis; analytic hierarchy process 
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Research on Heavy Metal Pollution of Topsoil 
WANG Jiu-ya, Ш Da-chao, ZHU Huai-jin 


(School of Mathematics & Computational Science, Sun Yat-sen University, Guangzhou 510275) 


Abstract: An effective method for locating the pollutant source and a categorization of pollutant sources are 
proposed. We first obtain spread pattern of heavy metals in topsoil through analysis of real data and simulation 
results, With different initial values, wo discover two different types of pollutant sources and the corresponding. 
spread patterns. One is constant pollutant source and another is abandoned pollutant source: Then we locate 
‘eal pollutant source by comparing the angle between gradients of floor surface and density surface. Finally, we 
categorize these pollutant sources by using surface fitting with Gaussian distribution and exponential distribution. 
‘based on the least square technique. 

Keywords: heavy metal pollution; angle between gradients; diffusion model; constant pollutant resource; 
abandoned pollutant resource 
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Contaminant Sources 
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Abstract: This paper presents a mathematical model to predict the location of the contaminant sources based 

‘on the profile of the pollution level of the against region. То do the prediction, we assume certain 

sampling sites to be the sources and simulate the pollution level of the region. Numerous math- 

‘ematical tools are used to simulate the diffusion of heavy metal through soil and air. By varying 

the assumption of the sources and comparing the simulated situation under different assumptions. 

with the real situation, we obtain the optimal number and location of the sources. Computer 

Implementation of the model gives a reasonable result. 


Keywords: contaminant sources; location; heavy metal pollution. 
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1 Introduction 


Nowadiys, heavy metal pollution becomes increasingly prevalent in our daily Ше. In some 
areas, it has n fairly negative impact on ош health. To reduce the pollution, we heed to 
accurately locate the contaminant sources. In this paper, we offer a pragmatic way to analyse 
the pollution data and predict the location of contaminant sources 

14 Outline of our paper 

1) Section 2: concept of pollution index; 

2) Section 3: theoretical fündament of our model; 

3) Section 4 and 5: results, conclusion, assessment and generalisation of our method. 

1.2 General assumptions 

1) The density of heavy metal of the sample site is the average density of the square 
kilometre; 

2) The land shape can be approximated by the smooth surface connecting all the sample 
sites; 

3) Only diffasion through soil and diffusion through air are considered in our model. Other 
types of propagation may exist, but are not statistically significant; 

4) The source of contaminant is located at one of the sample sites. 


2. Pollution index 


‘To measure the level of heavy metal pollution for each site, we adopt the Nemerow pollution 
index? 

Single contamination index Р, = C,/S,, where P; is the single contamination index of heavy 
metal pollutant i, C; is its real concentration, and 5; is its regional background value, 
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Nemerow pollution index?! 
P. ет, 


where P, is the Nemerow pollution index of the site, Pmax is the highest single contamination 
index among different heavy metals; and Pyye is the average single contamination index. 


3 Mathematical model to predict the location of contaminants 


In this section, we present the theoretical work of our model. 

We will first assume one of the sites in the region to be the contaminant source and simulate 
the pollution level of the entire region under this assumption. Normal distribution is used to 
approximate the distribution of heavy metals in a flat region. Several correction factors, such 
ав the altitude difference, the smoothness of land path, the direction and strength of wind, are 
taken into consideration. Then, we compare the simulated pollution level of the region with the 
real situation provided in the excel file to test the degree of similarity (using mean square error), 
We then change our assumption about the contaminant sources and до the simulation, After 
all possible assumptions are simulated; the assumption with the greatest degree of similarity is 
concluded to be the location of the contaminant sources. 

3.1 Sub-model one: diffusion through soil 

As the rea] world operates, normal distribution is the most general objective distribution 
law. ‘Thus, it is reasonable to use normal density function to approximate the distribution of 
heavy metal in a flat region 


1 
В) = s P (< в. 


where A is the assumed contaminant source, В is the point whose level of pollution is to be 
simulated under the assumption that A is the contaminant source, dA, is the distance between 
А and В, and /л(В) denotes the pollution concentration at В. 

Since the region under investigation is far from a flat region, we need to take into account 
the relief of the terrain between A and B. Thus, we introduce two modifications to the above 
formula: land smoothness and difference in altitude. 

Since the rougher the land path between А and В, the more difficult for heavy metal to 
propagate, we use the correlation coefficient r to measure the smoothness of the land path and 
modify the formula as follows 


ЕА(В) 


т) x falda). 


‘The correlation coefficient is calculated by using the point B, the contaminant source point 
А and all other site points whose z, y coordinates lie near (<1 km) the line connecting А and 
в: 


Sc Wah) 


үдо- а-а, 
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Неге we assume there're n site points lying near the line AB. Then we get the projection of 
these points on the vertical plane containing А and В, and denote them ав 21, гонг ‚т. For 
zü = 1,2,- ,n), hy is its altitude, and d, is the horizontal distance between this point and 
А. 


Besides the above modification, we also need to consider the difference in altitudes of А and 
В. Here, we need to use slop k and the oblique angle 0: 
ha-ha 
dan 


, @ = arctan(k). 


Since we just compare the final difference between the simulated pollution and real pollution 
to predict the location of the contaminant source, we can just use an assumed function that 
has the same properties of the real distribution function to calculate. 

As far as it's concerned, when the terrain between A and В is a straight line and А is above 
B (B < 0), the steeper the slope is, the wider the pollutants are able to spread. Accordingly, the 
larger 0 is, the greater ор will be. Besides, when (8 < 0), ов > о. The analysis of situation 
when А is lower than B (0 > 0) is similar. What's more, when (0 = 0), ов = о according to 
the definition; while when @ goes to infinity, the pollution will not spread upward, and thus 
|. Then we can modify ов by 0: 


ов 
ар 00088, 820, өз-- , 0<0. 
в | 820, ов Я 


Here о is the standard deviation of the normal distribution where the land is plain. 

The ов function we define obeys all those properties mentioned above. So it is reasonable 
and well-defined. 

‘Therefore, we obtain the final distribution function of the heavy metal pollution propagated 


Ма вой: 
а 


РАВ) = х Se ew 


3.2 Sub-model two: diffusion through air 

In this sub-model, we still assume that the distribution of pollution level can be approxi- 
mated by a normal distribution centered at the contaminant source. However, we need to take 
into consideration some other influencing factors: wind strength and wind direction, which may 
change the shape of the distribution. 

First, assume that the direction of wind is completely random, and then the distributions of 
the concentrations of heavy metals will follow normal distribution, which centres at the original 
point (contaminant source), and thus z ~ М (0,6). 

‘Then based on the analysis of the spatial distribution of pollution provided in the excel 
files, we find out that the distribution of heavy metals is severely skewed towards the left- 
bottom comer of the graph. Thus, it is reasonable for us to assume that the direction of 
wind follows the direction of the vector (-1, —1). We assume that the strength of wind is 
random over а long period. Thus, the wind will shift the normal distribution curve. Hence 
the mean of the final normal distribution is approximately the uniformly distribution: и ~ 
U (0, Laverage strength); Where Leverage strength is the distance that the wind with average strength 
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can blow the substances containing heavy metals forward in the direction of the wind. Thus, 
we сап calculate the distribution along the direction (—1, —1) by stochastic analysis 


d 
teln) = vm GE. 
Then 
Y (ЕЧ 
хә зше 
where. А 
вие) = уж, «ёй. 


Then we need to know how it will change the mean and variance of the normal distribution 
mean and variance which show how the wind inffuences the location and spread of the pollution. 
We attain this айп in the fallming way: 


E(X li o) = (NÛ = i 8) u = po) = 


Hence 1, 
E(X) = Е(Е(Х |) = Е(Е(М(и, и) = E) = "eene, 
Then 
Var(X) = E(Var(X |) + Var(E(X | и)), 
Var(X |u = по) = Уа(М(и = no, É) и = po) = 8°. 
"Thus we have 


Var(X) = E(Var(X |) + Ма БО |ију = Е(62) + Маки) = 63 + — 


The above analysis is based on two points which lie exactly in the same direction as the 
wind vector, namely e. Then we will generalize the aforementioned case to any two points. 
Suppose that the contaminant point is А(54, ya), another point is B(zy, js), and then there are 
two vectors, namely: e = (-1,-1); АВ = (zi — ra, № — №). 

‘We can then calculate the angle between these two vectors: 

„(АВ-е) 
Саве), 

Knowing the angle between these two vectors, we can obtain а new random variable which 
suits the general cases by changing the uniform distribution ji ^U (0, cos à Laverage тенен) and 
ме plug it into the previous equations in the following way: 


БЕІН Peep | + ЕЧ 1) вух) = Сенче co a; 


Тима strength 008 ОУ 167 2 


сова 


F(a) 


Var(X) = E(Var(X | u)) + Ма (БО) = 266) + Уа) = #2 + — e Ола. 


ОО 
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ву using this method, we сап caleulat the concentration ofthe ending point by knowing 
the contaminant point given that it diffuses only through air 

33 Combination of the sub-models 

From the: above two sub-models, we obtain the distribution of the beavy metal pollution 
propagated via вой and the distribution of the pollution diffused via air. Now we combine these 
two distributions to get the final distribution of the pollution. 

То caleulate the final distribution, we assume a fixed proportion т of the pollution is prop- 
эмей via soil. Accordingly, 1 —7 of the pollution в diffused by sir 

‘Assuming А is the source, we want to get the concentration of the pollution at B. Firstly, we. 
neû to obtain the ratio of the pollution concentration at В to the concentration at A diffused 
Ма soil and air respectively by the two sub-model. We use Олт, Саз, Ci» Свз to denote the 
Pollution concentrations diffused via soil and air at A and В, respectively. 

‘According to the sub-models, we have 


p 


Then, the final ratio of the total pollution concentration at B given А is the contaminant source 
ку е 
Са 


Hence, the estimated pollution concentration at B is Са = х Сл, where CA js the real 
pollution concentration we collect at А. 

Tn ош model, we use Nemerow pollution index to represent the pollution. Hence, the formula. 
will be Pf = p Pa, where Pj is the estimated Nemerow pollution index at B, given A is the 
ontaminant source, and Ра is the real Nemerow pollution index at А. 

Ву implementing the above method, we can get one dataset of the estimated Nemerow 
pollution indexes өг pollution under the assumption that а ‘certain point is the contaminant 
Е нее. Forinstance, we assume point i is contaminant source, and then we get Pj which 
Жена the estimated Nemerow pollution index of pollution at point Е given í is the source. 
So we can get 319 such Pis. 

"To compare the estimated value with the real one, we use the following formula 


E 
E (Pi - Pa)? 
= 
39 7 
where Py is the real pollution concentration at point k. 
Applying М same method to every poit (habging we mentioned in the last paragraph), 
жс сап get 319 datasets of the estimated values, and we also get 319 MSP: 


MSE, = 
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If there's only one contaminant source, we shall find out min, (MSE;) when i =l; 
point { is the contaminant source. ^ 
If there're two contaminant sources, we use Pj (i, k = 1,2,:-- ,319) generated above 
calculate the sum pollution concentration caused by two contaminant sources. For 
же assume points 1, 2 are two contaminant sources, and then we get the estimated: 
concentration at point k: 
Py? = Ph + Р]. 


Similarly, we get 


MEE‘ 319 4 
and find out miny {MSE; y} when i = m, j = n, and points m and n are the most posible 
two contaminant sources. 

By using the same method we could get the most possible k contaminant sources (k = 
1,2,3, ). But we cannot increase k to 319, which will be too large for programming. So we 
just increase k until the minimum MSE of k — 1 assumed sources is smaller than the minimum 
MSE of k assumed sources. Thus, we get the most possible number of the contaminant sources 
and the location of them. 

By combining the two sub-models, we can predict the location of the contaminant sources. 


4 Results & interpreting the results 


4. Results 

According to ‘Evaluation Soil Contamination’ published by U.S. department of interior, it 
is reasonable to assume that тд, = 0.38, то = 0.18, ro, = 0.64, то, = 042, пи = 0.21, т = 
017, пњ = 0.08, zz = 0.35; and we assume that о = 500, 6 = 5000. The results of the 
estimated contaminant sources and the related minimum MSEs are listed in Tables 1, 2, 3. 


‘Table 1: Four possible contaminant sources 


Four contaminant sources Minimum MSE 
(8209, 6018) (2383, 3617) _ (13694.2357) (15248, 9106) | 1892324 
(2383, 3617) _ (15248, 9106) (13094,2357) (13797, 9621) | 2496879 
(15248,9106). (13094, 2357) (9819, 6799) (2383,3617) | — 2754345 
(13694, 2357) (13797, 9621) (2383, 9617) (8290, 6018) | — 3964291 


Table 2: Five possible contaminant sources 


Five contaminant sources Minimum MSE 
(2383, 3617) (15248, 9106) (13694, 2357) (2708, 2205) (3299. 6018) | 1094.358 
(2383, 3617) (15248, 9106) (13694, 2357) (2708, 2295) (13797, 9621)| 1732398 
(2383, 3617) (15248, 9106) (13694, 2357) (2708, 2295) (9319, 6799) | 1921.301 
(2383, 3617) (15248, 9106) (13694, 2357) (3290, 0018) (13797, 9621)| 2457.625 
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Table 3: Six possible contaminant sources 


Six contaminant sources Minimum МЗЕ 
(2383, 3617) (15248, 9106) (13694, 2357) (2708, 2295) (3299, 6018) (13797,0621) 1108231 
(2383, 3617) (15248, 9106) (13694, 2357) (2708, 2295) (3209, 6018) (9319, 6799) | 1321387 
(2383, 3617) (15248, 9106) (13694, 2357) (2708, 2295) (3299, 6018) (8629, 12086) 1431216 
(2383, 3617) (15248, 9106) (13694, 2357) (2708, 2295) (13797,9621) (9319, 6790) | 1442418 


42 Interpreting the result 

From the mean square error, it is seen that when there are five contaminant sources, we 
achieve the minimum mean square errors and their locations are (2383, 3617), (15248, 9106), 
(13604, 2357), (2708, 2295), (3299, 6018). This is reasonable because all of these five points 
are frequently shown in the above tables, indicating that they are of a large probability to be 
the contaminant points. 

However, since the mean square error is still very large, we have to suspect our result: 
Because we have made an assumption which states that the only possible contaminant points 
are those points shown in the data file, but in reality, it may not be correct. If we revise it a bit 
and try the middle points between (2383, 3617) and (2708, 2295) and also substitute the middle 
point of their heavy metals’ concentrations, we can obtain a better mean square error which 
is 190.134, 5 times smaller than the previous mean square error, when we try 3 contaminant 
points. ‘These three contaminant points are (15248, 9106), (13694, 2357), (2545.5, 2956). 


5 Additional information & evolution model 


If time factor can be taken into consideration, we could complete our work to develop 
an evolution model for the geological environment of the city. Hence we need to collect the 
concentrations of heavy metals at different times at the same points as we choose previously. 
Analyzing the dynamic data, we may obtain the accurate о. As time goes by, we are also 
able to collect the weather data weather of the city, so as to analyze the impact on pollution 
distribution by rainfalls, the wind direction and the information about the rivers across the 
city. 


the information mentioned above, we could modify our model by correcting the 
pollution extent and wind direction in diffusion by air sub-model, adding the factor of the 
rivers and rainfalls. Furthermore, we can develop a new model to study the evolution models 
of geological environment of the city on the basis of the current model. In the new model, 
‘we could analyze the pollution distribution of а certain point as time goes on, instead of the 
pollution distribution in а large area at a certain time. This new model will also help us predict 
the trend of the pollution distribution as time goes on. What's more, we could predict what the 
environment we reside will be in the future, and the significance of the environmental protection 
will become a highlight of our life. Therefore, we are able to avoid the trend of the environment 
deterioration, and make our environment better and better. 


88 CHINESE JOURNAL OF ENGINEERING MATHEMATICS. VOL. 28 


and then 


т‏ و 
point 1 is the contaminant source.‏ 

ТЕ there're two contaminant sources, we use PL (6, k 1,2,::- 319) generated above to 
calculate the sum pollution concentration caused Му two contaminant sources. For example, 
we assume points 1, 2 are two contaminant оиго; and then we get the estimated pollution. 
concentration at point k: 


Similarly, we get 
эз 


AE - Pî 
ма. 
eid Sid йм BEng) when ть = n and pointe т and n are the most possible 
two contaminant sources. 

By using the same method we could get the mos. possible k contaminant sources ( = 
1,2,3,7-). But we cannot increase k to 319, which will be too large for programming. 80 we 
just increase k until the minimum MSE of k — 1 assumed sources is smaller than the minimum 
MSE of k assumed sources. ‘Thus, we get the most possible number of the contaminant sources 
and the location of them. 

By combining. the two: sub-models, we can predict the location of the contaminant sources. 


4 Results & interpreting the results 


ал Results 

According to “Evaluation Soil ‘Contamination’ published by U.S. department of interior"), it 
is reasonable to assume that za, = 038, 704 = 018. ros = 064, тс = 042, пу = 021, ты = 
0.17, ть = 0.08, па = 0.35; and we вешто ‘hat о = 500, 8 = 5000. The results of the 
ышана contamiñant sources and the related minimum MSEs are listed in Tables 1, 2, 3. 


ble 1: Four possible contaminant sources 


Four contaminant sources Minimum MSE 
Taago, 6018) (2989, 9617) _ (13604, 2557) (sa 9106) | 192324 
(Qm зт) _ (25248, 9106) ~ (13604.2857) (13190 9621) | _ 2498879 
"азма, 9106) (13094, эз5) (6819, 6700) (06% зв) | тыз 
desir) (19707, ет) (2389, 3817 (649. vou) | 3564251 


Table 2; Five possible contaminant sources 


Five contaminant sources. ‘Minimum MSE 
asas, 3617) (15248. 9106) (13094, 2357) (2706. 2295 (3209, 6018) | 1094258 
(азва, 3617) (15248, 9100) (13004, 2567) (2108. 2295) (13707, 9621) | 1732398 
(383, 3617) (15248, 9106) (13004, 2357) (2708. 2295) (9319, 6799) | 1921301 
(3383, 3617) (15248, 9109) (13094, 2387) (8898. 6018) (13797, 9621)| 2487628 
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‘Table 3: Six possible contaminant sources 


Six contaminant sources [Minimum MSE 
- (2383, 3617) (15248, 9106) (13694, 2357) (2708, 2295) (3299, 6018) (13797, 9621)| 1108.231 
| (биз, 3617) (15248,9106) (13694, 2357) (2708, 2295) (5299. 0018) (9819, 0799) 1324387 
1, 3617) (15248, 9106) (13694, 2357) (2708, 2295) (3299, 6018) (8629, 12086)| 1431216 
- (2983, 3617) (15248, 9106) (13694, 2357) (2708, 2295) (13797, 9621) (9319, 6799)| 1442418 


42 Interpreting the result 

From the mean square error, it is seen that when: there are five contaminant sources, we 
‘achieve the minimum mean square errors and their locations are (2383, 3617), (15248, 9106), 
(13604, 2357), (2708, 2295), (3299, 6018). This is reasonable because all of these five points 
(ме frequently shown in the above tables, indicating that they are of a large probability to be 
the contaminant points. 

- However, since the mean square error is still very large, we have to suspect our result, 
Because we have made an assumption which states that the only possible contaminant points 
° those points shown in the data file, but in reality, it may not be correct. If we revise it a bit 


nt of their heavy metals’ concentrations, we can obtain a better mean square error which 
190.134, 5 times smaller than the previous mean square error, when we try 3 contaminant 
nts. These three contaminant points are (15248, 9106), (13694, 2357), (2545.5, 2956). 


Additional information & evolution model 


If time factor can be taken into consideration, we could complete our work to develop 
evolution model for the geological environment of the city. Hence we need to collect the 
ntrations of heavy metals at different times at the same points as we choose previously. 
ng the dynamic data, we may obtain the accurate о. As time goes by, wo are also 
"ble to collect the weather data weather of the city, so as to analyze the impact on pollution 
distribution by rainfalls, the wind direction and the information about the rivers across the 


Given the information mentioned above, we could modify our model by correcting the 
ition extent and wind direction in diffusion by air sub-model, adding the factor of the 
and rainfalls. Furthermore, we can develop a new model to study the evolution models 
Of geological environment of the city on the basis of the current model. In the new model, 

could analyze the pollution distribution of a certain point as time goes on, instead of the 
ution distribution in a large area at a certain time. This new model will also help us predict 
the trend of the pollution distribution as time goes on. What's more, we could predict what the 

onment we reside vill be in the future, and the significance of the environmental protection 

become a highlight of our life. Therefore, we are able to avoid the trend of the environment 
iderioration, and make our environment better and better. 
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Research on Heavy Metal Pollution in Urban Topsoil 


ZHENG Ji-yong, ZHAI Dong-ju, JIN Mu-hong 
Advisor: LEI Hong 


(Department of Mathematics and Computer Science, Tongren University, Tongren, Guizhou 554300) 


Abstract: Aiming at solving the problem to analyze the intrinsic characteristics of heavy metal pollution in 
the urban topsoil, this paper studies the spatial distribution of the element concentration, the main causes 
underlying heavy metal pollution and the propagating character of the pollutions. First, based on the sample 
data, we obtain the concentration distribution of each element by using Kriging interpolation method. And the 
spatial distribution map is further drawn by virtue of the Surfer software. With reference to the background 
values of concentration, the degree exceeding the stipulated standards of the various elements ік attained in all 
functional zone based оп the average actual concentration by relative contamination rate. In addition, we assess 
the environmental quality with respect to the geo-accumulation index and provide the pollution levels of each 
heavy metal element in different functional zone. Moreover, the distribution difference of various elementa is 
‘analyzed based on the concentration distribution map. Based on such environmental difference, it is evident that 
human activities are the main cause of heavy metal pollution. Second, we establish the convection model for 
describing the pollutant spreading by virtue af the propagating character of heavy metal pollutant. АП possible 
locations of polluting sources are calculated reversely using the characteristic curves of space coordinates. After 
analyzing the intrinsic principal underlying the model and evaluating the experimental results, we put forward 
some further research orientations, such as the increase of the sample size, the improvement of the interpolation 
algorithm, and the establishment of the time series evolution model provide some new ideas, 

Keywords: heavy metal pollution; the index of Geo-accumulation; Kriging interpolation; convection model 
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зеңі чин. $, CSR ERR SRE 99 


СУЗЕ | $8 SFA k W о Ë ES sP A ТЕК HA (г.и), BRON by 
Airam) НА, B Z 0, ВАНЯ, PEERI 


dag = V (#1 — 22)" n — va x 0.1 (km). а) 


2) ЖІЛІ Floyd Wik) RAH EPSA REEI РОВ. 
Жа С БОЧНЕ Ao, НИЯ, RETI SR T f; WAIL 3min A fl 


EH O А И 28, 29, 38, 39, 61, 92357 ИН. 
412 RAR 
) ОРОН. fy WEERA ТАТА РЈУ УК 
асау (2) 
2) ранио. Чокаи LOU РЕЧЕ, № 
20 
Уљ=1. (3) 


ч 
3 BURT POR. UT Y Gm. мн TEN, ЖЖМ 
АНИ, ш 
fa =1. (4) 
4) FRAME 3min PREM AED AL, ПТ ОЛУРА SE, JUR СЕНЕ ЗА 3min s ff 
ЖОШ, UAE 9808098809 3km, РОА 


Јо x dy < 3. (5) 
5) МЕЛ 3min ARIRNI, ВАЖНЫЕ ОИС ее ЖИТ. В) 
do x fo < min(dy). (6) 


4.1.3 Hur 
жн EOE A AG ERO RON US 


LU очен m 
2021 


= 
w= M (fo xn), 
Е 
=, 
2 
Жизь 
at f faci, ® 
fo xdg <3, j# 28,20,38,30,61,92, 
fu x diy < min(dij), j = 28,29, 38, 39,61,92, 
fa € (0.1), 


1,2, „92. 


100 в & ъв ost 


414 HEX 

RUT ICL Matlab ЖҮНІ Floyd ПЕН, ЖӨЕ ЖСН AZ 0080 PROS DES MU 
RII Lingo (ACE НЕ T fr ЖИН. НОКИИ Rh ВН 
Spr бн Жж y А 3.02. 

42 жи. онна Лнан ЫШ 

1) MUS ar 

апаванавпн, AS UG SUUS GEL RI А, По RI тт 
dM es ЛОВИ MEF SU UE OS M. (ій, ЕЖЕ, ЖАН 
Бийн иа kat EF HF BHF ТАИ РОВ 


LEE Lg 1 
min max ( fy x“). (9) 
воли 
v а, 
ша 3 3 fy x => (10) 
=“ 
ПЖ: 
m 
Жіті. 
E 
да fe sy (11) 
fo € 40,1), 
i= 1,2... ,20, j= 12,14, 16,21,22, 23, 24, 28, 29, 30, 38, 48, 62, 
2) Кн 


ЖИИ, РЕЈ EBS Ж EM Е 8.02min AS СДАНА УС WAT Н А РА В ВОИ 
MN. 

аз WI. MAMET EPROM 

1) МАНИ 

Б ЕИ RA HIE 3min Py + AE (ІҢ X T 3km) ЕГІЗІ 28, 29, 38, 39, 61, 
925648 O. EF HEE 9128, 29 9 ë O REFEN Рав МП 28 9 #129 ВО 
по да 9138,39 9 Pš L10998 Fë Т T 3km (R8 1141 38, 39, 40 FFE: ДЕ ВК 61 $ 
пй (BREAST % F 3km ПЯ 48,61 9 ЖП: #2502 9 ИЕП AN 
F 3km Г 87, 88,89, 90,91,92 92871. БИЖ, EHE GRE HERE ИХ: 3min ІН 
хай, CEREAL A SBE BSA GERA в, REME НТ УА З РИО 
иж, тиан. 

2) МЕХ 

EAU B RII ЛК, ЧЕТ а R m KF RST МИС, 4 
ПЕШ ЕГЕР. 

3 кажат: 


min mx (хт). (12) 
ES 


Li тан. % РАИ SAIR 101. 


AIR 
m 
тш 3S (o.- 9). (as) 
ена. 
fj € (01), 
di x 1,53, 
а 
т. (14) 
3) Масим 


T ADIRI Lingo ESR 9 — ЕЛИН 
=. 
max (У f wj) = 85. 
m 


БОРО RI СТО РОЗЕ ви У MP, PU Lingo СЕЈ ВЕРЕ ОН 
JY GER. CF AZ ВЕ аа 29 99 Пу, A22 V f: fo W GE 
WAGEN, A23 3668393648 ОП, A24 ВОИ 87 ИО. 

НАР SY f: t Sñ E E ЖЖ 00773609220. 5 E WW G2 E 
003,02 НИЕ mei, ТЕШЕ mA F a. 

4) узнава 

ШЕШИНЕ ERATES DUE HH: А Es УЖ ЖКА та тв 
AALS FTE OCF, (АЗР а РАА РЛЫ Н Е, КО 8.5, КАМИНИ 
ЖИ 5 НИТ k. 

FEDERE, UPL, RM, ПЕНА нито G A25 AREN 22 PCT, АНЯ 
JP SOUS OT Co EI EIA ЈУ E 1.51, MK THEMEN 6.1, МА T HON 4 NE 
frt 2208.5, ИЯ P Виан жия. 

4.4 Мад. PHRA LERRA FAEERE 

4.41 Я У FT OOM AMF 

1) МАРА ХНК AK, WHAE, HDR MPRA AMER, 

2) CRAMER PHF ARMM, ВИТ, 

3) ОБЖ ТЕЛАХ, ADS, Пеене, 

4) 412 T AK 4F 6 3min f ELE UA OB RR > ICH a BE У G 
НОВИНАР 3km MORE, BRAC, CRAG, DEKH 12^, EKHI, ЕК 
извг, ЕУР ЈЕО СРЕСКЕ CRDI AMR. 

442 би ЖЕШ та 

1) get 

UTAH BUR HERRE тт ласына КОШ. AT as HR ва 
Af Ж ЕКИП SE, BEARS RE. 


102 ов м 4 2 4 E: 


ЖИП ЖН БИ E BAY ЕЕ RE FG. ROUTE JE UC 
f ЖЕЉЕ НО DU ШЖ, (UI БОК RUS АЖ: АХ 0.62634, ВК 0.38877, CK 
Ж 0.842281, DKW 0.502818, EK % 0.807227, F I % 0.69527. 

ЖЕЗ ОКИ ИЛЕП A Us A PC 15.9%, BIG 9.8%, C P£ 21.3%, DK 15%, EK 
204%, F EG 17.6%: ERY 6 8600980, РОВА ИНК К ИНЖ 13,8,17, 
12, 16, 14. 

2) PARER a v 

CT 8E BC XE CY. f; A W (> ЖИЙ E DL ДР P fo А IB DO IZ пи 
WI0-1 WERBEN AHT RL REIL. ОПАТЫ В С E I АНИ. 


нй... 
питке 
min max ( s их ш). (15) 
E 
ERI: 
Ju € (0.1), 
kee (0,1), 
m 
Қ-8, 
ча 
СЕ е EM 
ж Ляо Mk=, 
E 
мек 
3) атдюжи 


АПИЯ Lingo ЖЫН, RAY SW YG CACO. ша Bi. 

45 MEENET 

1) Миних 

РАПН ле AERE AGAR IU Ас ЖИН, РЕЧЕ IESUS ЕРИ ИСА 
MU E IUE А FTE NE свит ср Н ЭШИК ЖИН, AEDU ls ВИ ЗИ A E 
Ж. запетая РЖ, Меги ЕУ E HUS EA ЕЕ ЙН A ТА О 


S) В ay nin 
өйи Fei ARE, Е И “- 


3min. 
Жл EERE, И 
ПАН, Зе: 7,8, 9,30, 31, 32, 33, 34, 35, 36, 39,45, 46,47, 48,61. ЯНА, RIII 
ЛЕВА, В (А, B #@2 11066867): 


А = (ААА). PEA, _ В=(ВуВ;-- „Веј, q€B, (18) 


тнт "ЕЕ. % RRBs FS OR SECURE 103 


Jib A ARE AFR, Rate ERE PARTE TROLS. ОТЕ ВО AU 
НОО НЕ: BAe BITE. Reem AM OZ а, of 
ЕКОО OR AA Me SRA. А A, B UG e T RR Ө k, 


Жа: BALHPRESEFERAR 


STARE ЕГЕТ жеее бюз к (BORG) 
А 13 1,3, 5, 10, 11, 12, 13, 14, 29, 31, 40, 60, 91 
в 8 94, 95, 97, 08, 100, 103, 125, 152 
5 x 167, 168, 169, 171, 174, 175, 176, 177, 178, 
180, 181, 200, 210, 239, 262, 275, 300. 
+ 5 307, 328, 331, 332, 339, 343, 344, 360, 362, 
363, 368, 370 
š $ 372, 373, 374, 377, 379, 381, 383, 384, 385, 
387, 388, 390, 396, 404, 419, 457 
1 " 47, 478, 479, 480, 482, 483, 484, 485, 505, 
мї, 549, 574, 579, 582 
2) ВВ 
Не. 
(19) 
жт. 
(20) 
#=1,2, 
3) Нижа 


ЖЖМ, ОНЕ АВ 22715 НИИ, account (A) = 45, IPRA dig № 8.79 
km, ОИЕ ñ: 8.79min A SEE E P MERE А EMH, ЕВУ AO E EEE 
тап. 


5 Ше 


ИТЧА ГВА, MAA АЕ ЗВЕНЕ ЕЖА 

1) ERE ВИАН Р, RIUTAN ҖЕ} F€ ЖУ бокпуъ, A ЕП 
4, шина ЕУ ЖЫ: 

2) МИЛАНО ЖААТ AH, HE EMM, ГЕК НИНИ 
ERK. 


104 r . . oW 22 


важи 


П PD. поема усетим Lk: SWT НЕ. 2000 
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[a] Жән, WK. Matlab ЖШ ЯГ SEIM]. 40» ВИД ЧИН, 2009 
Song Y Z, Jia D Y. Matlab Numerical Analysis and Application|M] Beijing: Mechanical Industry Press, 
2009 


Research on Setting and Scheduling Scheme of Traffic and Patrol 
Police Service Platform. 


YE Yunjia, LIU Jian, WANG Yu 
Advisor. САО Hua-lin 


(Aeronautical Department of Mechanical Engineering, Naval Aeronautical 
Engineering Institute, Qingdao Campus, Qingdao 266041) ; 


Abstract: The paper investigates the distribution situation of the trafic and patrol police platform in a certain. 
city, We develop а reasonable programming for trafic and patrol police service platform that needs to raise 
the resource utilization rate and enforce Ия response speed by using the Floyd algorithm in grapy theory as 
well us the least square and multiple objective programming. Furthermore, we investigate the unreasonable 
factors underlying the current police scheduling condition, and by combining with the practical situation we 
further propose the advanced scheme that re-echedulés the police strength. Besides, we establish the scheme of 
blockading crossing roads based on specific situations. 

Keywords: Floyd algorithm; least square method: multiple objective programming. 
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X44 8:1005-3085(2011)07-0105-12. 


IRERE FARE SSD 


н в, њеве, вит 
(MAREMARE, UK 401311; 2 WEFA THER ETHER, MDM 450002) 
а З жооп кеш Енен в сако жавокова” на. AAT 
MAESA, кпанентихледвниожалиндунаки, DHE E RW ION 
ж, AMT Geo ATED RAR. RGR TOME paymi Hm. 
XU. козина, man manu тщш 
##8: AMS(2000) 9030 MARS: 0221.2 дана: A 


1 Те ска У 


«ЖШ TR 9 F f; D ЕЕЕ АЕ Ж Ж ЖИН 
аве, + FJ poti аса ЕО GM TITRE. 

и кзоожо тн, ЖИЛИН, ЮНИТ 300 4 283€ 
TOF A, 8300 £087). МА "ЕЖЕН. HEME. хапни, ВЕЗИ" ИХ. 
ЖОНИ КК. f LA IEE FHK. TORN “ЛТ”. QA 
Өкен Е, аве, GIS. СРОК. НИЩО. BAVE iH 
Ж. Зо. OEI. WR UW. EMH, O 4 AW ATO. 
PETIG ВЕРНЫ “хин” жа. СРИ А ВУ ИЕ A 
ERREI, эшн Ап, BEEMAN. РЕМИ, Wakin. НО. 
ж. ROMA ВИ, кекаюжи. жанас, MT ATE, EEM. 

REUS DUE жи T 110 122 fü) fr REN, Qa Y T ПАКАО ЕНИ). ВЕЩ 
RORUR A it E (0 5 A 10 y hiii А 3 305 y bk. {ЕКЕМ таа, НЕЗ RAE 
HS GRE. WEWEEE COBRA, KURM AH, возия. Я-А, Ж 

МЕ”. 

жалға ЖАН НЕ ЛУ ВРХ. MERRER, ВЕНЫ, ЛЕН]. EMET 
# LEA dE. BAMA M TRE. Жат жн. Vif ga NOSE 
du. КЕ “Хит” жж, ккиняжаана. WP 28 ид. ща. FUER 
д, WDR AAE, WAPARE. НЫНЕ галон, HR COH HAE. 
пана. кана, ШЕШ ОЕГАН ЕК Ж. H TEA e SM 
тв, WEE) анине о ESOS EUN y. ТЕЛЕН 
. inf fth ORAS S OR GRO ER GE JR FA кт IDEST. | 
Hr И TREE. 

ЖА IE RITE 3E 2010 F 4 H E ЖПК ЕШР ЕИ ОЕШ Я, чөй A 
жы. възнавтихо- ева. иш: пш К (f fni erae, rft Hn 
а, BA CARE (FA женен? — MERA ES Dib? (OE AZ ЕНИ 
ат, ХВЕ? RA ae лгі? AG. 742010 8 AAA 


106 т а ж 4 ан most 


джан, RAREST AANT, ISLET ZS SSD М, CORT кат 
ЭИ ЕРНСТ НО РАСЕ РЕ ЧР. 11342010 EMITTER PE: 
JA LT MMR. жЕ ОУ MRR, HEE ORTE 
DMT, {12011 #5 AER ТЖЕ И. — КЕКЕК 
[ср Жр ыша у ИВА. 43 


2 ARR MARS FA KE 


ЖШ А ВА ВЕНА FS RU 

Татан PETE, MERE он a is 
deti CU 

2) FOIE аз ORLA HR. ER ERIS Y Л); K 

3) НИЗ, GE ЕЕН взе е RCA КИИ ИЛ; Ж. 

24 ABFA MMT M A AA REM. 

211 AR SCR RERO FOOSE ЕҮТ? 

MEA CIA HH Хоа ESMER. a Ж ЕД = L2, m) ИИ, уу 
FMC = 1,2, n) BERHE бт. пипат. ФО АИ, КАМ 
ЭЛИ НИ LA A M — A ОМА ИЕ ШИН ek 
HL = (нна BBARR AE ЖАНЕ ВОО Floyd ЗАВ, Л Matlab Я lil Y ЖК 
W Wa ZONE, GIS NE А, (ij = 1,2. ут), MYL kn. 

2.1.2 АОИ 

ГТ КЕЙТН. WIL А Т а ЯНЕ Та 
fikir. RMU REPE ER ЧЕЙИ DZ ii: 

1) ЧЕЗ RE НЕГИГИ 
ТАНДЫ а 

2) диня, жететін. 

АЕ, ОЧИТЕ Х = (Хун. ЗЕРЕН 

1, FE yy TN; 
"ode xe 
Hl Bí ЖЖ АМИ z, BF у, (i = 1.2,--. m. 22 п) ЖАН (ICM fW) 
= (Ву). ФОЛИ X F. OCA AAE LS TUI Y e AD ЗЬ РЕ (ЧЕ) E: 
их 


= (Ту)=х» = (Xu ` Bü)mxn- 
та, ЕР ЖИ, ЖАН 


зах Ту, jolie. 


По ТИ АМИН ЖЕНЕ AHE W = (uiia, лет), Hw; ARAKA 
RO (G= 1,2, m) RRR MEF APRA AAS = W -X = (61,63: Gn): 


Lun нң. 3. REEF TOMER ARE 107 


BT CEPA LRT, AES RATHER, WER 


a= |> (6,-8Y 
| $ 
Mih, КОРЕН. 


бела, АКИТА К fr EAK EATER AEREE В TAE 
НЫ, ibo deor e GI A В RNF: 


min mue mino(G), 


Ex i-a 
ЈЕ 

Ж.Ха>і, jelAhegn 

а 

Xgy20Xl HÜi-jM, Xy=1, 
d-212..,m, )ж1,2,-.-,п. 


S25 Ae PIC PE EA CAE EL OLE ОИНИ, З fk А 
‘Sh CORO, WETAN, НОМ D il тая. 

213 PAM SER AT 

T LTA AL Вж O-1 мит. PEOR ILA. Зак Е, РОА 
UH, ATO AE {ЕЗ АЕ ЗЧ). Ш 28-45 29,38,39,61 81929081, За, НОЖ. 
ШИЕ ТЕР 159. 168, 28. 7920 FER. KIEZ A O REET CL A 
зн Ж, WREKE 


st. 


max Ту $3 
За“ 
Амир, SUNOS АД FE EF тәрік (ЕЕ) JI 78 FL PE 0-1 ЗЕМЯ. 

5) —JH o WN H REITER ATR A. ЧЕЙН T 61 ЗА ЖЕНИХ 
йк. ROTH ЕНИ ВДВ ЖИН ЕТ. РЕЧИ СИТА РОДИТИ 
ГТК Н КК КІРЕ. M Lingo КӨДЕ A. UB Matlab FER WW ЖГ. 
SUR dq. мотив, MBRR, WREE RR. 

REPT y ay RHI ЩЕ: 20 A ft P ES ACCU R Wy 2.1759 PH, ВАН 5.7 y 
bb; JAPU ЯН 6.225 c. НЕНА TSK, MDM LG NK: ЕЕЗ АИ АНА 
dieti. 

22 AROS S THEBRDIS HAS. 

M 20 437-6 rt 13 TOREM URS H bak CURT — UBER UR (dij aos RR А 
WEG л (i = 1,2, , 20) AERO z, (у = 1,2,--- ім (НМ) ERE. FE yı RE 
FO 2, ATHURI ЖЕНЕ X = (eyo: ЖУАН у HERO SM. zo =1, Ж 
Ети = 0. 

ЧЕХ, АРЕНЫ ОВ (НЮ) РЕМ T = (Гоа = 
(Фулу)лоз, ЗЕЛИНА А РУР) max, Т). kB tA Н ЛЕШ КИЕ Л, юн 


108 r в а t * а _ mou 


шинж DASARI ке: у 


1 
sg =O, Ван-1, 
i=1,2% 20, 3 1,2,- „13. 


Зеро вот аж А вп, КЕТ БАЊИ. 
яве тен. Jl Lingo РУСИ СЕНА, RYH LI EI Ж ЖИП 
МӨНІН, SLE 


Ki лая (HALA: #8) 


те Tia Тани е [Haso [нане | z€ | ияле | ниче 
кеа» pres [na [asss aspa aa Рае [s 
м [| оз [me| | ws [АШ] ж | чи 
lw [xem an] м [sme me| м pe 
~ | э | sos |a| 12 kag š 

a | № 1 зове |Аз| 2 | 05 


изаа 980155 944. 3H HOM 3.55296 FH. жәе, Xx 
ARD IT KK REKE E -in. 

23 ARMIN EMERAN 

«аи, кино ебати emma, нажил» 
EER: HENEK HF Er ESATO RAN. PO ЕЕ А 
KH нин T {ЕНЕ БЕСЕ, ГҮҮ wR HB. BAR MMT 
м HASH ви. 

EMA FERRE, ООВ Amn MEEK т — n PT 
ГІЛІ EAE лай) КОЛЕ, авит ева V Gr DAE ULT Ж 
Loi м-н P GOD de tnnt ях. dcr se an Y f 
MERIR. HA Es psu ол, мината. =4. RA 
звания иа. 

REET BLOF В OE BET HABA ха, тан ^^ ат Wanu ЛИ 
позвана А (rns ra ST) зри ERE Gro 1r Xen = 0. 

Agittin ЛЕВ. ПОЗАДИНИ RORIS FAA п +1, BB) MATER, XE 17 
жашан ERAS, заря ОКНЕ зада, ЖБИ id 
LLL OR UL 6 2 Ni fi Y EEA SERE to A. AE RAI ЕЛ оо. eM 
жакша свежа) A 

tot 


к=" +) 


99—01 
ж” 


Lun вн. % заветине 109 


жөнеле жн. та. TEFAF E KRR: 


со 


b-t 
пах Рр +(1-9) 


т =08 1, i=1,% m-n 


ЭМЗ, SIENE (1) ОЕ, ОГО АЗУ AS ORR ERR. 
И ORG, поето fr. RERE. Ж = 05, НЕМ 
ILE SAS G ТАИ ТА R. Ine 2 BOT. 


Жа MIRE GATOR (HARA: e) 


жекен. | изн или | хижа | нажена | unm | нжх | лежа. 
хжале | тю | 2000 | - | RSF | 3682 | 24157 | 0045802 
нате | 4102 | 29545 | 0077239 | Mika ке | 36822 | 19062 | 0086828 
иа ње | 36822 | 26503 | 0110575 | Mestre | 34071 | 1.8426 | 0.048671 


отт, ТЕЛЛЕ GI AGAS c тйк bitte 
Лу АЖ НИЈЕ 29, 40, 48 #187 ОП. 


з SHRM MS FAM 


ЮЛАЙ ЕК ОШ A F iñ 4 48; 

1 WORDT A BEY RM RHE, питижая, WEHM R: 

2) KIL PRAET ООЖ, BERE Hh UIE TF fF VERD ЛИНЕА Ж. 

за треш жао 

cl I у АУ, ЕЕ A PR WG — 
AEA IIH: TTA — о ла ын И 
ЛЕННЯ. Д и HC, WPW IE 1) HIER T 

XP AIDE GR аа A тан ВНУ 


T= max (6). 


- n) ЖЕ В y, FRE BA (GL RC IUS IHE FG. ИО EDS GG = 


в-а ол, 


А5060) 


WJ TP 6 ана 122222 


по шн & 4 9 а ase 


ники TER i TERRA A DC AE, PL RR ER 
м, BIST CLAS P 6 ИИНЕНИН ЕО Ж. Апия, EANN 5в2 НО o 
BIM 138 ЛЕЛЕ ЕЗ PARE, ЖЕЙН ван PRESSE а. HITE RG OT E 
SIM IURIS ново ft F SERIO УЖ, FEREN 80 FE RE 
ERI RI 3.7735 949, ЖЕСИН ЗНН ПЯ 12.6803 FH: ЖЕРЫ ЖИН Ж 8.4312 ЖК, 
ЖЕНЕВА, RDA LOK. вися, ЭП аласа 
qma eto. 

32 шоителжнажен. 

AE RP OUR TEE ERA. LAA BY RH ER 
шуй. магнити A KOREK. АОИ, ASME ATR CURT 
PARKI, RATE Cr вт EKE DETER f. Вай, AER ML XC 
ЭНН ИССА ВИННИ, WA М AA A TOME HE 
ЗРЕНИЕ, НОННА ST ОВС. ЯРАР z ДЖИ Е 

t а, 
т жаа жа бы 
кер 
t= Уи = = Жағу 
т 


ЖИ т, НОСНЕ CES 


ОТГ НЧ А G, WWW т, жетті, 


Gmax = max (С) 
BUR e ЖСН AD BURG IEAA E ӨЙНЕ RIQ (h, n 
жакт КМГ ЕМЕ “PEAR, НВ EE ИАТА W ATB NI A ЈЕ 
Mis {ЕК BR VA EGER ЖЕКЕ ЖЫН. как взвод, а ЕН E GAP fr tt 
[MILLS 


Аз SUG Hota oT 


митек | mite Gis | акжени (24) | RATE (RIA) | HEHE 
° ЕЛ 126803 261 44946 
1 за 12.6803 184 41365 
2 517 126503 184 3.9251 
з зав 9.8992 184 3.9753 
4 330 8.1188 184 39973 
5 288 8.1188 159 зата 
6 206 8.1185 159 Л 


тиз эя. з. жажнятвнакзиянт ш 


Жазаыян» MST FONE ERS REE E ШИСИ Ж, НЕЙ 6 4-5 
НЕ В ИКЕ. AE НН 5 РЕА СРЕ Раса RU 
Жак. 

зз ФЛЕШ ЖЕНА ERS 

ЖИАНЫН Р ARE Bs ПР 5 ROS ЙЕР (ЯНИ) ЖЕРЕ D = (Ds, Da, --- , Dn) ЖЕ. 
ЖИВ ЖЕШ ВИДО ORBEA i$ Pa 00 2F RE EH) бо km/h. E REET, ЛЕ ЛОЖНО 
ЖИ, TREE AYR RAL ED ЖИ. ЗОНА Ној АНЕ DLN Р Do rp co RC IL 
T) PARURE САД. ZEND t PER арс ML VA PTF BE В 21004625402) F(t), HOS MORE ЛО 
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Optimized Model of Setting and Scheduling of Traffic and Patrol 
Police Service Platforms 
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Abstract: In this paper, according to the problem B of 2011 Higher Education Press Cup CUMCM, the 
background of the problem of setting and scheduling of traffic and patrol police service platforms is introduced. 
‘Then the approaches and the concrete models of several related problems are given, as well as the main solving, 
methods and the corresponding results, and the general existing problems of the participating papers are analyzed 
briefty, Finally, several problems worthy of further study are given. 

Keywords: traffic and patrol police service platforms; the transportation network map; model of distribution; 
surrounds model 


Жов ана вз. ‘Vol. 28 Supp.1 
2011F11 | CHINESE JOURNAL OF ENGINEERING MATHEMATICS Мон. 2011 


38 5:1005-3085(2011)07-0117-08 


Еа FARES HE 


MLS, жең, жез 
поеми: FRY, KTH, Rie 
(жт PIER, W 713800) 


жан. ишин, ‚ ктт. пищевом али. ELIO SOS. «иш 
ИНФ, КАША XUR, RU Matlab 1 Lingo КАНОН ДИНА W. FEZ A Ret =f IN 
ЖИВ — iie. китет. KARRA HEYA N. 

я ШАРОН ЬО ЕИЗ Р т. езана. A E 
TIO — AUN. AGRUM Floyd Wik. WE A ЧО AEAEE fr. KE f 
{AVOID АА КО A НЕ хит. АННО, сини 
SIDES ан, ИТЕЛЕ Ча ARD H NAR, тиетш. (ULTIO 
ж, N 138 AR ER PPA Se cs, BLAS 9, 1E DR E fr U AIR 
AUS MRE 6 PEK ОИ УСНИ ИВЕ 
RM BUR ASA вита 2, 


жий: олш, Floyd Wks OBE, mE 
94%: AMS(2000) 9029; 90820 AMAR: ош! змей A 


1 вази 


11 Б-я 

НЕ О БЛАНКА ТА СОЕ ЕЈ РАТЕ А а 
(АЕ ЛОИ ЕТ, JE DAT A не MGR CIE сини O ОНИ AGE OE FAN, AHIR 
ЗЕ ВЕНЕ а VIERGE ACE OR. 

2) ОКИ ЕАН ІЗ PARR ME, AEE ИН RANK K 
Медија НИН SOIRS IM Ж Gr BERI ТОДОР. RE EG PAE ER UC 
db, UE ETRE н SF HEEE URS 

3) no EM WERE 3 y OA А SE RDU ERIE EMRE, РАСА AL ЖЫ 
H, Уве AE ФЕР, Redi HERE CLR (АЖ) мята: завзет HERO IU. 
BIHE T fE R) 7736 Eo i ROSE RERO. ЗВЕР АЕ GO C Ж. 

12 NA 

1) AERA AEE HF Gs IBN REA. заявления OV 
ВИШЕ ЖИВИНЕ. DUE. ELEN IER HUE REX IURE IH AT IH DR 
E. {ЕЗТН Ж Жан E PERO GEI ӨН a Ж. 

2) втажян, птшилесват, AIR аи. ик, «#онт 
м, жалық “ннен PCIE AF НИЕ А UR DL IO РЕН 
IM 3 8” ФЕД НЕВИНИ ПА ACIE. НЕ ТЕК ЧА ННН. 


из таж та зове 
2 впюво5жи 


FRA, МАНЕРЕ = 60 km /h ӨН LTE и”, ВН 
кийа 


тою = !0тт/ш, 
НН ФР PPPH RIN Е 10mm (РА. ЖСҚ, ANEH Floyd MN Matlab 
Mat REH FESP fa S Йй 

24 ERNIE SRW 

211 GO ION 

RR НА N RAER, ЖИНА БОНН TC ARCA СРЈ 9 P f; t: 
ж. E RF REOR AP OE A ly, IER GROS 


ГА = mindy, j=1,2,.-. 


HEA KEEFE ИКЕМ МЕЛ; Ж, Е PR: 


AL ADREPFCEMARM FEAR 


же | жала +6 тияк [te] чили 
1,67, 68, 09, 71, 73, 
' 2 | 2,39, 40,43, 48, 70,72 | з | 3,54, 55, 65, 66 
та, 15, 76,78 
5, 49, 0, 51, 52,55, 56, 
а | 4,57, 00, 62,63, 64 АУМА АШ 
58, 50 
т |е 92, 47, 48,61 | в | 633,46 9 [931,34 35, 45 
10 [1 п | 11, 26,27 HEE 
з [uana | а [n 15 | 15, 28, 29 
16 | 16,97, 38 їт [ата 18 | 18, 0, 8I, 82, ва 
SERE за | 20 84.85.86, 87,88, 89, 
90, 91, 92 


RY AE, LL Ht RIL A KHR 2925, 29, 38, 39, 61, 9296 IGUR E 
TEE 3 PHA ІА. BEK ROE OE EY РАДИ 6 PEL КИЙИН, (nde 2 Bo. кн 
ЭД, Ө EIERE CEN F é E: 


Жо REF EHR GAR ABO HRH tu: ##) 


же [se | annm | же | њо | анны | же | зо | анима 
16 | ss | 3408 | 7 | е | 41e | 2 | з | 36829 
2 | 92 | 36027 | 15 | 28 |-алам | 15 | 29 | 570063 


Lun жт. 9. завет AN 19 


212 аннан нити нат жи 

1) B THEBSE ST ми 

ағашы, ПНЕ ЖЖ Fs ЗУ, (20 GERRE РАН, 
OM 13 ЖЕНЕ ЗЕМУНА, LAAN. NLM, Жо КР, 

a) EO  XERCREI ЖР GI 131, 13 46238 НИМ BEL. fq P| ТИ, 

b) нат жым 13 вина, EUH IMR T JUPE A 

с) иажаквинашнтен, МЕНЕ НИМ TH ft MER. 

tide e еве 


min e= mm Бе) 


Mah 
„йана j 120 XU, 
0, hi TOF CARN VOR SOON, 
ly RAM ОШ REP В ВВ j OR AE АЕ. с) = aul, ЖЕНОМ j ЖОШ ЖОН 
инини. 
BIRR A 


20 
Pagel, 3-12-25, OR MACH ZA dI ЕЖЕН ЖИ 


= 


14 
Уау, i= 1,2,0 ЖИРНЫЕ НИИ HT 
я 
UH Lingo РИЙ е = 80.15 mm. 18 GË 9 R Ра f t FT ZE 8.015 2 OWN 24 13 Ae 
ШОТ еме, LEE WI RM RIME. ЖЖ, ER H Е A ЊЕ 
Jy 461.88mm. 


аз делата 
7523 тања E ялан 
2-3% 2-ю—з0—3в 4— 62 1-6 
5-4 = 7-7 1530229 
830 | 8—33—32—7—30 | 9-16 | 9—35—36— 16 
10—22 10——11—22 r= 1-25-24 
12-12 | Rer 1323 КЕРІ 
uon мз 15—28 15—28 
16-14 r= 


120 I “Ж ТЕ . жй gast 


2.1.3 MMF ARS (THRRET) MEL SRR 

HEBD ARGH Juk Tr НЕЕ ИЕ EF E ЕЈ ЖЕР, By «30mm: 

EEK НИЖЕ T fF ARE FY ЕН ЗА ЧЕ fu] T RIER SCR ROS RL) T” 2 ножа 
Ж, WAAR ET ERD W, = P (T, + Тур. ПЕЙКИ E QU, В 
A ЕТ 


АЕ. аиту, 
emp FN 


Ie 2 
min 52 = = > (W, -WY 


іші 
Wis X Тур, 
530, 
2<п<5. 


HEDL Е ВЕНЕ P ЕУ SU 6 2.1.1 I II GA AS A У ОВИСИ ТУ 
Ff ИЗА ОГ IMN, HERRE HET REIR GU LR. DEOS £ 
ЖИМ (28,29), {38,39,40}, (48,61), {87, 88,89, 9D, 91,92) аскета + 
WELUN GLA IE Д КИ. MUR СЕ ЗЕН, за Matlab женена, OBI ALDER 
HOP fi ft: W 28, 39, 48, STIA. ЖЕН JC HERE КИ, 05.6185. I cC i L 
TE UL EE КЕЛЕТИН, WOE А Tr fW UTI. Жо ЛЕЛЕ 
Mh t 001890 V fa t W 69. ОК гана талу жал ТОГА M T ИУ 
3069142925. МОИ. Dt, BIRE ERREKEN, ПОТИ 
WS ЛЖИ OP fr. АЕ 28, 39, 48, 87 XI Wi C IA REN UO НЫ Во ЗАКОНИ i) 
WH. 

Ht fi BS Dt Е ЛИ P r feia a ВИ, Inde 4 ЭГ. 


Жо REG SEM 


Fe] жала Fe] Fuam #6 гала 
1 TOOT „лшщ з | з |3,54 8605, 66 

74,75, 76, 78 
4 | 4,87, 60,62, 3,64 | лды |a 

56,58, 59 

т | в [Геза ° | 9,31, 34,95, 45 
10 [16 п | 11,26, 27 HEE 
1з | 13,21, 203,24 HH 15 (a5, 28,29 
16 | 16, 36,37 irae 38 | 18,80, в, 82,88 
19 [197,79 20. | 20,85, 86 28 [25.28 


39 | 38,39 as | 47, 48,61 87 | ва, 87, 88, 89, 90, 91, 92 


ттт мош. з. TROP GU RR SAA 121 


22 AM RENET SRE 

221 пятезианнен 

ДОРТ Е И Е, H Floyd ERAT A PFE ЗП RO R I КИЙ 
пан. EXC EE 138 AOR EX EH ИЕ ETCH BI, PMA ИЖ 
fHE REX ЕЖ. ЖЗ PP fi Е У мж Ж. аа 
молнкак, ERAF E RRR лама (7; Ж—) Жалел ежен 
MEYE (HR) PAG AER A BE fs 

ШУ T it tg 4" Ратови (РИТА P T fr СО ОСТАО R kf N 
дати. 

222 тирани 

MT 22.1 ғ. ИМЕНИ Г ЖЗ: ЖАН ЗАРЕ ТАЖ 
м. 

а®— EARLY Ee Mi MIN TOR UY ет IT ALB СУЗА ЕЛ 
ik, Уре ТКЕН. Я 


min а-ә 
= 

4 $30, i=1,2,... (а) 
Ens ј=12 (5) 
гу ср, #12, © 
zg-0Hl, i=1,2, m, (4) 
ю=0@1, #=1,2,т, (e) 
рті, #=1,2, 0 


Зета DETER, ЖЕЖ ВЕЗАНЕ, Бъгита ща 


га 1 ома, 
0, жтвожвев. 
планов, ly Ri FOI j EO РО 


1 Dee EM у “ип, 
у= 
47 o жое тап. 


(a) Str EASY ВЮ 5k 0 Pr fi а езана p|: (b) os SE R Н AP fF 
в. (с) RANA TER ТИП УШК GIA BERE TROLS, (D) Ren ЛИ 
пене. 

ката ГРО, ЧЕН Рае BUR OF ТЕ 
dS. FH 5.1.3 PTH НОН Ж. 

Fi Lingo MERGE FAR. MK SPUR. 


122 


тщш mis. = ж uan 


жу MARR 


ga | tret | tea 


25, 39, 48; 87 


104, 149 


184, 199, 201, 204, 205, 216, 238, 240, 252, 257, 259, 263, 388, 304, 318 


329, 332, 333, 335, 344, 362, 369, 371 


387, 389, 390, 392, 393, 408, 419, 420, 421, 446, 454, 458, 472, 474 


»|s|v|o|w|» 


491, 509, 513, 517, 539, 541, 558, 560, 569, 574, 575, 578, 582 


MLZ жаве, NRT AT, WN I о ЕЕЗ FP r Oe 
BH, ЭМЕК КЕ T fF MM, ROM, вее 


nis S 


= 
путу 530, i=1,2 m, P= Lym, 
Dayal, је от, 
ача іші2,-,т, ДРИН m, 
zy=0R1, i-1l2--,m 
р=0 1, +=1,2, m. 


Ml Lingo @ AML Y AHR, 1056 РЖ. 


Же мажат 
EM | умљељи | инея 
А 15 3, 10, 11, 14, 15, 22, 25, 20, 33, 40, 43, 48, 58, 66, 90 
в 5 104, 114, 135, 149, 163. 
н = 166, 175, 177, 179, 184, 185, 193, 199, 201, 204, 205, 214, 222, 231, 
238, 240, 250, 254, 256, 259, 263, 266, 274, 285, 303, 310, 319 
D 12 327, 329, 332, 333, 338, 344, 350, 361, 362, 363, 369, 371. 
5 22 373, 378, 387, 389, 390, 393, 397, 400, 408, 416, 419, 420, 421, 423, 
424, 440, 444, 446, 457, 459, 470, 472, 473 
x ч 479, 485, 488, 489, 499, 509, 511, 515, 532, 539, 541, 559, 562, 565, 
570, 573, 574, 575, 582 
наше Бато, йж давюотжа, пат ПЕЛЕ 56. ША 


13646. AW ЕВО ЕНЕЛ, EARN С ЛЕНЕ ЖЕ. 


ана моз. ¥: тшшен кая 123 


2.23 пзвааназивс5жи 
жишш 2.1.2 PREJAKA, PALES итн, ту 
RIN M80088 Ps UI ТИНЕ AGERE ajak 3 UU 


(еее Los (Ye) e 


Vet RT ЙН! F: 
(ТИЈ 
. mint, 
PURE 
2 
{min e = | шах (55 еи) <-%, 
vt < rg, 
= 
ag =O 1, 
м 
Lay =1, 
я 
Lay <1. 
я 
UH Lingo WERI ALA CEI Ж. nde T Bos. EL Kh, RIMKE 7.96min ІК 
арғын, абе, HEAR. 
Ат RR, 
AS [muse [43| mana AS | sui 
1 [122,4] 2 [атат 3 | 12, 27, 26 
а | 15, 28, 29 5 вота в [i225 
то | 88088 |y | обв 9 | 10, 34, 35, 45, 44, 2, 70 
60, 62 
10 | 19, 77, 76 1r |а, 62, 190, 180, 168 12 | 173, 232, 211, 216, 215 
із | лт, 208, 218 | 14 | 172, 227 15 | 169, 254, 253, 240 
16 | 167, 249, 248 | z | Po 224 228. 225, 18 | 320, 349, 371 
170, 273 
jg |i 491, вв, | 2 | 47 ав, 495.483. 122, | „| | 476.545, 846, Ыт, 532 
487, 489, 482, 491, 486, 487 548, 549 
aa | 9595555 | щ | ers 557, 564,563, 562 
557, 558 


124 T m oe ж а gusti 


3 Hum 


1) rU ESSET. EA TEFLA Matlab, Lingo Ф, BMT DEAF. 
2) ивента АНАЖИЖ, AA RRM. 
чивтатетедвашищо, DD CIAM 


Las d 


li] хиж. ЗЕЛЕНЕ AIM). Жй: PRT, 2003 
Wu Q Z. Operations Research and the Optimization Method|M]. Beijing: Mechanical Industry Press, 2003 
[2] Rafik, WW. Matlab Жі to S [M]. EH: магия, 2009 
Song Y Z, Jia D Y. Numerical Analysis and Application of Matlab[M]. Beijing: Mechanical Industry Press, 
2006. 
ls] MA. MM, көт. Matlab 0 аяйм]. k, МТН, 2008 
Lin X 8, Zhou Q, Lin D X. Matlab7.0 Application[M]. Beijing: Mechanical Industry Press, 2006 
№] жим, шем, не И (м). Ло, азаннян. 2004 
Jiang Q Y, Xie J X, Ye J. Mathematical МодеЙМ). Beijing: Higher Education Press, 2004 
lr] Wik. ati ni RPM). ERO. ирип. 2005 
Han Z С. Mathematical Modeling and Application|M]. Beijing: Higher Education Press, 2005 


Setting and Scheduling of Traffic and Patrol Police Service Platforms 


HU Lizhong, ZHANG Jin-lei, ZHU Jin-bao 
Advisor: LI An-ping, AN Qini, LIANG Fang-chi 


(Missile Institute, Air Force Engineering University, Xi'an 713800) 


Abstract: It is an important issue to allocate and schedule the traffic and patrol police service platforms 
(TPPSPs) in accordance with the practical police force and condition. By employing the Floyd algorithm, the 
distances between the locales and TPPSPs are obtained abd the jurisdiction region for each TPPSP in the city 
center District A is allocated, Based on the object function of minimizing the blockade time and the shortcut, 
the blockade scheme for the TPPSPs of District A is scheduled. The supplemental TPPSPs' number is accessed 
‘on condition that the police can arrive at the locale within 3 minutes with the minimum workload variance 
of the TPPSPs of District A and the minimum number of supplemental TPPSPs. Through investigating the 
reasonableness of the TPPSP setting for the entire city, the 138 acces arteries cannot arrive at the locale within 
3 minutes, The models of added and reallocated TPPSPs are adjusted on condition that the police can arrive 
tat the locale within 3 minutes. An optimal containment strategy is established according to the object function 
of minimum blockade time, region area and region population. 

Keywords: 0-1 program; Floyd algorithm; short path problem; balance degree 
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‘The Setting and Scheduling Model of Traffic and Patrol Police 
Platform Based on 0-1 Programming 


TANG Dil, DONG Xisoxiao', WEI Xin! 
Advisor: LUO Wan-chun* 


(1- The 9th Team of Cadet Brigade, Third Military Medical University, Chongqing 400038; 
2 Department of Mathematics, Third Military Medical University, Chongqing 400038) 


Abstract: This paper aims at resolving the problem of the setting and scheduling of the traffic and patrol 
police platform. First, the shortest path and distance between the intersection and platform is calculated based 
‘on the Floyd algorithm, and each intersection is set to be dominated by the nearest platform. In emergency 
‘ese, the model is established on the basis of 0-1 programming to attain the fastest and the shortest schedul 
scheme. By compromising the workload of each platform and the time wasted on the way to the intersection into 
consideration, we establish a partition model for the territorial jurisdiction and achieve the scheme of vetting 
new platforms. Finally, in order to catch the criminal fled from Point P as soon as possible, we construct the 
optimization model by setting the likely longest feing time to be shortest as the objective and that the police can 
*ucoesnfully close the intersections and the criminal can be completely blocked up as the constraint conditions 
‘The best plan for blockade can then be drawn by solving tbe model using 0-1 programming. The final result is 
ne follows: the police from 18 platforms should be scheduled to close 18 intersections and the criminal is to be 
‘caught within 20.52 minutes. 

Keywords: shortest path; 0-1 integer programming: the police platform 
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ЗК map, WHAM ЖҮ, Wi 5080096 CU E HC 28% MRS. 

PROBA MARAT, N tHE ILA nila fr АЈ нож, каним Ж ERR 
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FDU T HEROS FUR. WERNA 0р, ТЕЛЕРИ ЈА ЛО D SENI É 7) 
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ЭН mF RE TOS. JF AREER ORT: БІН, R F GRE 
DAR. 
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аз, РИА IE GU X EERE, RRD, KATHA: AEM. 
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Setting and Scheduling of Traffic and Patrol Police 
Service Platform 


XING Zhong:yang, ZENG Fan-zheng, ZHANG Jiang-bin 


(College of Opto-Eiectronic Science and Engineering, National University 
of Defense Technology, Changsha 410073) 


Abstract: In this paper, we investigate the scheduling and setting of the traffic and patrol police service 
platforms, aiming to optimize the patrol time and balance the workload for each platform. The Floyd algo- 
rithm is used to attain the shortest route, and then the workload of each platform is optimally compromised. 
Accordingly, each of their optimal administration area can be adaptively allocated and tbe location of added. 
platforms can be generated, that is, 28, 36. 48, 57. 81. Besides, for quick blocking of the 13 traffic arteries, 
the Hungarian algorithm is used to calculate the shortest time for a complete blocking, namely, 8.015 minutes, 
As for the adjustment of the platforms, we propose the evaluation function F, including the expected workload, 
the variance, and the longest patrol time. We finally evaluate the optimal scheme for besieging suspects within 
the blockade circle. 

Keywords: shortest path; cover within 3 minute; traffic and patrol police service platform 
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Modeling the Setting and Scheduling of Traffic and Patrol Police 
Service Platforms 


WANG Wei, II Xin, ZHANG Jian-peng 
Advisor: LUO Wen-chang. 


(Faculty of Science, Ningbo University, Ningbo 315211) 


Abstract: In this paper, we consider the problem of setting and scheduling of trafic and patrol police service 
platforms. Fist, we obtain the shortest path matrix of nodes by programming with the Matlab software, and 
then transform it into the shortest time matrix of nodes. For the frst question of Problem one, we pre-asign 6 
nodes that the polion cannot arrive at within 3 minutes and then propose the corresponding 0-1 programming 
models with the goals to minimize the maximal workload of police and the variance of workload of police under 
the constraint that the responding time of police is no more than 3 minutes, respectively. By comparing these 
results, we show that the latter is better on compromising the workload of police, For the second question of 
Problem one, we present а 0-1 programming model with the goal to minimize the longest responding time of 
police subject to blockading all the 13 acces arteries. By solving the proposed model, we show that the longest 
blockading time ік 8.0155 minutes. For the third question of Problem one, we derive the corresponding 0-1 
programming models with the goals to minimize the maximal workload of police and the variance of workload 
of police under the constraint that the responding time of police is no more than 3 minutes. Then the results 
ме compared, For the first question of Problem two, we propose three indexes, ie., the responding time of 
police, the workload of police and the total responding time of police, to evaluate the reasonableness of the 
present schemes, respectively. According to the three indexes, we propose the corresponding 0-1 programming. 
models мер by step to evaluate the reasonableness of present schemes. Based on the three indexes, we give 
the evaluation results for the six districts (A, B, С, D, E, F). By resetting the platforms or adding the new 
Platforms, we give the adjustment schemes for the six districts (A, B, C, D, E, F), For the second question of 
Problem two, we first get the nodes of each level from Р node to arrive at by programming with Matlab. Then. 
by using the iden of gradually giving back to contain the suspects, we propose the containment strategy. ‘The 
suspects can be contained within 7382 minutes by using the new strategy. 

Keywords: setting and scheduling of platforms; 0-1 programming: responding time of police, worklond of police 
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Тһе Optimal Setting and Scheduling of the Traffic and Patrol Police 


ZHANG Хо, TIAN Zhi-xiang, FU Jian-long 
Advisor: Instructor Group. 


(School of Finance, Shanghai University of Finance and Economics, Shanghai 200433) 


Abstract: This paper first employs the Floyd algorithm to preliminarily allocate the police stations, and builds. 
the network-flow model to establish a scheme to blockade an area. Second, a multiobjective programming model 
is used to evaluate the current allocation of the police resource and present the corresponding improvement. 
Finally an encitcle-catch model is set up. The network-low model is utilized again to make different encircling 
plans according to the suspect's speed, and on the basis of the greedy algorithm the encirclement. under two 
constraints is constructed to catch the suspect within one hour. 

Keywords: network flow model; multi-objective programming model; simulated annealing algorithm; compre- 
hensive evaluation model 
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ЕВИЯ E Matlab РЖИ, SE F: 
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жин “аж” ва, EREE E ILIO HU A BE J PU НЯЛИАЫ ALIENI ЗЕЙ Ж. 
цента. RAGA ACRI CUR ERE НИНЕ Ж, ВЕРФИ. 
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clear; close all; cle 
zuo=xlsroad(’C:\Docunents and 

Settings \Administrator\/M И сиаси.х]в" , 1, 'В2:С583'); 
lvexlsread('C:NDocunents and 

Settings Administrator cumca.xls',2,'A2:B929'); 
dist=zeros(1,l);k=l; 


dist(1,k)=sqrt((zuo(1v(i,1),1)-zuo(1v(1,2),1))"2 
*(zuo(1v(4,1),2)-zuo(1v(1,2),2))"2); 


kokel; 
end 
istedist/10;1u=zeros (582,582) ;k=1; 
for 1=1:928 


1uQvG, 1) ,1v(4,2) mdist (k); 
1uQv(1,2) 1v(5,1))«distQO ; 
көкі; 


if заа, о 
лаб, 3)=10000; 
end 
end 
end 
for i=1:682 
лаб ,4)=05 
end 
weight=1u; 
for k=1:582 
for + 
for j=1:582 
if weight (i, j)>weight (i,k)+weight (К, j) 
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weight (i, j)=weight (i,k) weight (k, j) ; 
end 


tic speed=60; УМЕНИЕ (ка/һ) 
station=x1sread(’C:\Documents and 
Settings MdninistratorVitiliVcusca. x15" ,3, *82:8817); 
nunstation=length(etation); FHA AC 
spespeed/60; numA=0; 
for 1+1:582 Ж КОННА А ТАВА (РАН: 
if (weight (32,4) « Зевр)) 
numA=numA+1; 
А(ошвА) =; 
end 
ond 
taile1;times=0; 
‘while ((numA<583) ва (1а117=0)) 
times=tines+l; 
fa11*0;nunB*0 ;nunD-0; 
B=0;D=0;Dstation=0; 
station(:,2)=1; 
for 1-1:928 
df (Cásenpty (find (дтн (4,1), 1))BkisenptyCfind(Beelv(i,2), 1))) 
за (isompty (find (am1v(4,2), 1)))) 
nump=numB+1; 
B(nunB)elv(i,2); 
E 
Af(Cisempty(fina(A==1v(i,2), 1))&(isempty(find(Beelv(i,1), 1))) 
ttlisonpty(find(A==1v(i,1), 1999) 
numB=numB+1; 
B(nusB)-1v(1,1); 
end 


re 
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Ddistance=0; 
di-veight(32,B(1)) /sp; YD BIREW à f) RH [а] 
for j=l:numstation 
d2(j, D-veight(station(j,1),B(1)) 3 877 АВИВ ra И ib IN f 
ева 
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Def ind((d1-d2) .«station(: ,2)>3); 
if (isempty(D)) 
пшаА=пшай+1; A(numA)=B(i) ; 
B=[B(1:1-1),B(i+1:numB)] ; 
failei; break; 


numD=numD+1; 

for k=l:length(D) 
Ddistance(k)=42(D(k),1); 
end 
(Onin, nunDuin)emin(Ddistance); 
Dstation(numD)estation(D(numDmin),1); 
station(D(numDmin),2)=0; 
ond 

end 

end 

toc display(A) display(B) C=(A,B) 
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Setting and Scheduling of Traffic and Patrol Police Service Platforms 
LIU Jin, YANG Yun, WU Ke 


(Faculty of Mathematics Statistics, Wuhan University, Wuhan 430072) 


Abstract: Assuming that each police platform can only block one point, we create a dynamic optimization 
model for the optimal setting and scheduling problem of the traffic and patrol police service platforms under the 
consideration of the uncertainty of the criminal's Beeing direction. We first create a subset A which consists of 
the points with loss than 3-minute away from the criminals fleeing start point. Then we formulate а cut-set В 
‘which consists of the adjacent points of A, and search and allocate the adjacent police platforms to the points in 
cut-set B. If no police platforms can be allocated, we reformulate cut-set В again until an appropriate blocking 
‘condition tallies with the requirements. In this paper, we evaluate each nearly optimal blocking strategy when 
the criminal's speed is 60, 70, 80 or 90 km/h. Specifically, when the criminal’s speed is 60 km/h, the fastest 
blocking time is 8.9 min. 

Keywords: 0-1 programming; Floyd's algorithm: platform scheduling; best block way 
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Abstract: This paper discusses the problem of optimization in natural easing assembly. Based on the thorough 
study on its background and requirements, we propose two mathematical models and analyze their computing 
behavior in integer programming. 
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А Heuristic Method for Solving the Assignment Problem of Natural 
Casing Materials 


YANG Zhong, ZHANG Qiliang, XU Xiao-hui 
Advisor: LING Weiwei, Ll Fan-chun. 


(Department of Basic Courses, Jiangxi College of Applied Technology, Ganzhou 341000) 


‘Abstract: In this paper, we establish an optimization model for the assignment problem of the natural casing, 
‘and design a heuristic algorithm to solve the problem. For the requirement (1), we reduce the problem into 
a maximal number of bundles for size 1, size 2 and siae 3. Based on the Planning theory, we establish an 
‘optimization model for the problem on condition that the number and the length of each product are fixed and 
Known, and also we get the maximal number of bundle is 182. For the requirement (2), we obtain the material 
number in each bundle of the scheme, which gets the maximal number of bundles by employing Lingo 10.0. А» 
compared with all schemes, we obtain the best scheme which gets the maximal number of maximal length of 
shortest materials іп each bundle. For the requirements (3) and (4), frst of all, we obtain the number of each 
bundle size on the assumption of no degraded using, and then we get the bundle number of remainder materials 
on the assumption of degraded using, Finally, we get the maximal bundle number 208 and obtain the optimal 
scheme. For the time limit requirement (5), based on the establishment of large-scale integer programming 
model, this paper introduces a creative heuristic algorithm and obtains the approximate solution. The programs 
of the proposed heuristic algorithm can attain results within 1118 seconds. Worker can arrange the production. 
according to side ВИ в prescription. Also, we give an improved neighborhood point search algorithm to solve 
this question. We further compare the results and complexities of the two proposed methods, 

Keywords: casing assignment; optimization models; heuristic algorithm 
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Study оп the Maximum Number of Bundles about Collocation of 
Natural Casing 


ZHOU Ya, LIU Zhen-tang, LU Xue-yang 
Advisor: HUANG De 


(Guangxi City College, Guangxi 532100) 
Abstract: This paper attains all product schemes of constituting a bundle of natural casing from the company's 
specific requirements. Because the scale of the attained schemes is too big, we abandon those raw materials 
‘more than three kinds of above schemes, and build the bundle number maximum model in the presumption that 
the raw material cannot be relegated to use. In order to improve the utilization rate of the raw materials, this 
paper considers that the raw materials сап be degraded to use, and draw a general integer programming model 
for all the products. The optimal solution is calculated from all of the model solution bundles, The optimal 
solution is 191 bundles, and thus the optimum matching scheme can be determined. 
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